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^«d9e^Spplication of: Philip KUSK 



IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 



Application No.: 09/889,491 



Group Art Unit: To be assigned 



Filed: July 18,2001 



Examiner: To be assigned 



For: GENETIC PREDISPOSITION Attorney Docket No.: 8969-029-999 



RESPONSE TO NOTIFICATION TO COMPLY WITH REQUIREMENTS 
FOR PATENT APPLICATIONS CONTAINING NUCLEOTIDE SEQUENCE 
AND/OR AMINO ACID SEOUENCE DISCLOSURES 



Assistant Commissioner for Patents 
Washington, D.C. 20231 



In response to the Notification to Comply With Requirements For Patent 
Applications Containing Nucleotide Sequence And/or Amino Acid Sequence Disclosures 
(hereinafter the "Notification") mailed by the United States Patent and Trademark Office on 
August 31, 2001 in connection with the above-identified application, please enter the 
following amendments and remarks. Applicant submits herewith: (1) a Sequence Listing in 
paper and computer readable form pursuant to 37 C.F.R. §1 .821(c), (d) and (e), respectively; 
and (2) a return copy of the Notification. 



Sir: 



AMENDMENTS 



IN THE SPECIFICATION: 



Please amend the specification as follows: 
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Please replace the paragraph beginning at page 17, lines 10 through 23 with the 
following: 

Figure 1 panel A shows the location of the two polymorphisms, called BSP-A1496G 
(SEQ ID NO. 13) and BPS-G1869A (SEQ ID NO. 14), in the bone sialoprotein gene 
promoter. Panel B shows the location of the polymorphism, called MGP-C242A (SEQ ID 
NO. 15), in the matrix gla protein gene promoter. Panel C shows the location of the 
polymorphisms, called OPN-G520A (SEQ ID NO. 16) and OPN-T1825C (SEQ ID NO. 17), 
in the osteopontin gene promoter. The wild type sequences encompassing the four 
polymorphic sites for all three said promoters are shown with the nucleotide at the 
polymorphic position in bold and with the substituting nucleotide - also in bold - positioned 
above the polymorphic site. All nucleotide numbering is relative to base pair number 1, 
which is the most 5' nucleotide of each of the promoter sequences as published in the 
GenBank nucleotide database. 

Please replace the paragraphs bridging pages 22 through 24 beginning at page 22, line 
21 through page 24, line 20 with the following: 

DNA analyses. Screening for the BSP-A1496G (SEQ ID NO. 13) and the BSP-G1869A 
(SEQ ID NO. 14) polymorphisms (basepair numbering according to numbering of BSP 
promoter sequence submitted to GenBank, accession #L24756), the MGP-C242A (SEQ ID 
NO. 15) polymorphism (basepair numbering according to numbering of MGP promoter 
sequence submitted to GenBank, accession #M55270), as well as the OPN-G520A (SEQ ID 
NO. 16) and OPN-T1825C (SEQ ID NO. 17) polymorphisms (basepair numbering according 
to nimibering of osteopontin promoter sequence submitted to GenBank, accession #D14813) 
were performed as follows: 

The polymerase chain reaction (PGR) was used to amplify approximately 250 bp long 
DNA fragments of the BSP, MGP, and OPN promoters encompassing the BSP-A1496G 
(SEQ ID NO. 13), BSP-G1869A (SEQ ID NO. 14), MGP-C242A (SEQ ID NO. 15), OPN- 
G520A (SEQ ID NO. 16), and OPN-T1825C (SEQ ID NO. 17) polymorphic basepairs. PGR 
techniques are well known in the art and it would be within the ambit of a person of ordinary 
skill in this art to identify primers for amplifying a suitable section of the BSP, MGP and 
OPN genes including the positions 1496bp and 1869pb in the BSP promoter, the position 
242bp in the MGP promoter, and the positions 520bp and 1 825bp in the osteopontin 
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promoter. PGR techniques are described for example in patents US4683202 or 
EP0200362B1 . Two hundred ng of genomic DNA was added to 25 \xl reaction containing Ix 
Taq polymerase buffer with 12.5 mM MgCIjCPerkin Elmer), 5 nmol of each dNTP, 20 pmol 
of forward and reverse primer, and 1.25 units of AmpliTaq Gold (Perkin Elmer). The 
reaction was heated to 95 °C for 9 minutes followed by 35 cycles of 95 °C for 30 seconds, 
46°C (BSP-A1496G and BSP-G1869A polymorphisms) or 49°C (MGP-C242A 
polymorphism) or 46°C (OPN-G520A polymorphism) or 48°C (OPN-T1825C) for 30 
seconds and 72° C for 30 seconds - the latter incubation with a 5 second time extension per 
cycle. The reaction was finally incubated 7 minutes at 72°C for completion of the extension 
reaction. Primer sequences for PGR amplification of DNA fi-agments encompassing the BSP- 
A1496G, BSP-G1869A, MGP-G242A, and OPN-G520A and OPN-T1825G polymorphic 
basepairs were: 

BSP-1496G polymorphism primer set: 
Forward primer: 5' - GAA AAG ATA TAT ATA GAA GGG GAA G - 3' (SEQ ID No. 1) 
Reverse primer: 5' - TAA TAT GAT TTG ATG TTT GGT GGT G - 3' (SEQ ID No. 2) 

BSP-G1869A polymorphism primer set: 
Forward primer: 5' - TTG TTT GGA GAT AGT GAA AAG AGG T - 3' (SEQ ID No. 3) 
Reverse primer: 5' - GGT GGA TTG TGA GGA GAA AAG - 3' (SEQ ID No. 4) 

MGP-G242A polymorphism primer set: 
Forward primer: 5' - GAG TGA GAA AGG TGA TGA GTT GGT G - 3' (SEQ ID No. 5) 
Reverse primer: 5' - ATT GTG GGA TGG ATG GAT GGA TGG A - 3' (SEQ ID No. 6) 

OPN-G520A polymorphism primer set: 
Forward primer: 5' - GGG TGG AAT TAA GAA AAT TGG TAG A - 3' (SEQ ID No. 7) 
Reverse primer: 5' - GTT GTG AAT TTA GTG GAG GGA GAT G - 3' (SEQ ID No. 8) 

OPN-T1825G polymorphism primer set: 
Forward primer: 5' - GAG TAG TAA AGG AGA GAG GGG AGG T - 3' (SEQ ID No. 9) 
Reverse primer: 5' - GTA GGT TTT TGA TTT AGG GGA TGG G - 3' (SEQ ID No. 10) 
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Please replace the paragraph beginning at page 27, lines 15 through 27 with the 
following: 

Results 

Five previously unknown polymorphisms were identified by sequencing specific 
promoter regions from the human BSP gene promoter, the human MGP gene promoter, and 
the human OPN gene promoter following a PGR ampUfication of 40 DNA samples from 
healthy women. The BSP-A1496G (SEQ ID NO. 13), BSP-G1869A (SEQ ID NO. 14), 
MGP-C242A (SEQ ID NO. 15), OPN-G520A (SEQ ID NO. 16) and OPN-T1825C (SEQ ID 
NO. 17) polymorphisms were coded as Xx, Yy, Zz, Bb, and Ss, respectively, where the 
uppercase letter signifies presence of the wild type base pair at the given polymorphic 
position and the lowercase letter signifies presence of the base pair different from the wild 
type base pair at the given polymorphic position. 

Please replace the paragraph bridging pages 27 and 28 beginning at page 27, lines 30 
through page 28, line 9 with the following: 

Table 1 shows the genotype distribution of DNA samples firom the 18 years study for 
all identified polymorphic sites: BSP-A1496G (SEQ ID NO. 13), BSP-G1869A (SEQ ID 
NO. 14), MGP-C242A (SEQ ID NO. 15), OPN-G520A (SEQ ID NO. 16) and OPN-T1825C 
(SEQ ID NO. 17). The left panel shows the actual number of samples categorized into three 
genotypes for the 3 identified polymorphic sites. The right panel displays the same analysis 
as the left except that the nimibers represent the percent of total samples analyzed for each 
polymorphic site. 

wt = XX, YY, ZZ, BB or SS 

hz = Xx, Yy, Zz, Bb or Ss 

pm = XX, yy, zz, bb or ss 

Please replace the paragraph beginning at page 29, lines 3 through 20 with the 
following: 

The genotype distributions for the 5 polymorphisms are shown in table 1 . For the 
BSP-A1496G (SEQ ID NO. 13) polymorphism the homozygous wild type genotype was the 
most abundant, followed by the heterozygous and homozygous polymorphic genotypes, with 
the homozygous polymorphic groups being quite small. In the case of BSP-G1869A (SEQ 
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ID NO. 14) polymorphism, the wild type genotype, as defined by the BSP gene promoter 
sequence from GenBank, was rare, the heterozygous genotype was 10 times more frequent 
and the homozygous polymorphic genotype was twice as frequent as the heterozygous. For 
the MGP-C242 (SEQ ID NO. 15) polymorphism the heterozygous genotype was the most 
abundant followed by the homozygous wild type and homozygous polymorphic genotypes. 
In the case of the OPN-G520A (SEQ ID NO. 16) polymorphism, the homozygous 
polymorphic genotype was the most abundant followed by the heterozygous and the wild type 
homozygous genotypes. The genotype distribution of the OPN-T1825C (SEQ ID NO. 17) 
polymorphism was, generally, the same as for the MGP-C242A (SEQ ID NO. 15) 
polymorphism. 

Please replace the paragraphs beginning at page 32, lines 1 through page 34, Hne 16 
with the following: 

From the table it is clear that the BSP-A1496G (SEQ ID NO. 13) and BSP-Gl 869A 
(SEQ ID NO. 14) polymorphic sites, especially when combined, are good sites for predicting 
whether an individual is genetically predisposed for high or low BMC/BMD. The OPN- 
T1825C (SEQ ID NO. 17) polymorphism only has a marginal influence on BMC/BMD on its 
own. However, when the OPN-Tl 825C (SEQ ID NO. 17) polymorphism is combined with 
the BSP-Gl 869A (SEQ ID NO. 14) polymorphism the percent separation of genotypes is 
better than either polymorphism alone (Table 2). On the other hand, the MGP-C242A (SEQ 
ID NO. 15) and OPN-G520A (SEQ ID NO. 16) polymorphisms are, at first glance, not 
suitable sites for such a prediction. None of the identified polymorphisms appeared to have a 
statistically significant impact on the change in bone mass over time (data not shown). Age, 
height and weight of the individuals involved m the 18 years study did not differ significantly 
between any of the genotype groups (data not shown). 

These observations strongly indicate that the BSP polymorphisms influence peak bone 
mass rather than the rate of bone loss. To substantiate this an analysis of the variation of 
BMD as measured 4 times on the same individual from 1977 to 1995 for the different 
genotypes was performed. In 1977 the average age of the individuals included m the 18 years 
study was 51.1 years, thus ending at 69.1 years in 1995. The expected outcome of a plot of 
the means of BMD for one of the two BSP polymorphisms as a fimction of time would be 
two parallel curves, each representing BMDs measured m individuals with the wild type 



genotype and BMDs measured in individuals with the polymorphic phenotype. Figures 2 and 
3 show that this is, indeed, the case for the BSP-A1496G (SEQ ID NO. 13) and BSP-Gl 869A 
(SEQ ID NO. 14) polymorphic sites. Moreover, the two BSP promoter polymorphism act in 
concert on peak bone mass to augment the mean BMD difference between genotypes even 
more than the isolated contributions of each polymorphism (Figure 4). 

The role - if any - of the MGP-C242A (SEQ ID NO. 15) and OPN-G520A (SEQ ID 
NO. 16) polymorphisms in the MGP and OPN promoters on bone turnover was less clear 
from the first analyses compiled in table 2. However, when BMC values grouped according 
to genotype were plotted as a function of time a set of curves appeared suggesting that both 
the MGP-C242A (SEQ ID NO. 15) (Figure 5) and OPN-G520A (SEQ ID NO. 16) (Figure 6) 
polymorphic sites are determinants of rate of bone loss. It is especially noteworthy that the 
ZZ and Zz+zz curves as well as the BB+Bb and bb curves separate between 1979 and 1989, 
corresponding to an average age of 53.1 years and 63.1 years, indicative of a genetic 
phenomenon associated with the menopause. Like the BSP polymorphisms, the combined 
action of the MGP-C242A (SEQ ID NO. 15) and OPN-G520A (SEQ ID NO. 16) 
polymorphisms also leads to a bigger difference between genotypes than either would create 
alone (Figure 7). 

According to the results compiled in table 2, the impact of the OPN-T1825C (SEQ ID 
NO. 17) polymorphism on BMC/BMD was only visible after it was combined with the BSP- 
Gl 869 A (SEQ ID NO. 14) polymorphism. From a graph of BMC values grouped according 
to genotype and plotted as a function of time it is difficuh to tell whether this polymorphism 
has an impact on rate of bone loss or peak bone mass, due to the proximity of the curves 
(Figure 8). However, combining this polymorphism with the BSP-Gl 869A (SEQ ID NO. 14) 
polymorphism gave rise to a set of time course curves clearly showing that these two 
polymorphisms cooperate in an additive fashion, and, hence, that the OPN-T1825C (SEQ ID 
NO. 17) polymorphism may influence peak bone mass, as the BSP-Gl 869A (SEQ ID NO. 
14) polymorphism, rather than the rate of bone loss (Figure 9). 

Finally, the association between genotype and urinary osteocalcin (N-MID®), 
Osteometer Biotech A/S) as well as urinary collagen type 1 C-terminal crossUnks 
(CrossLaps® Osteometer Biotech A/S) were examined. As expected there was no significant 
difference between the mean value of either of the biochemical bone turnover markers for the 
BSP-A1496G (SEQ ID NO. 13) and BSP-Gl 869A (SEQ ID NO. 14) polymorphic sites (data 
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not shown). Also, no significant difference was observed between the mean value of the 
biochemical bone turnover markers for the MGP-C242A (SEQ ID NO. 15) and OPN-G520A 
(SEQ ID NO. 16) polymorphic sites. This is likely due to the equally paced bone loss of the 
ZZ and Zz+zz genotype groups as well as the BB+Bb and bb genotypes groups at the time of 
N-MID® and CrossLaps® measurement (1995) according to Figures 5 and 6. 

Please replace the paragraph beginning at page 37, lines 1 through 19 with the 
following: 

The impact of the identified polymorphic site on bone mass as represented by bone 
mineral content (BMC) and bone mineral density (BMD) measurements at the distal arm in 
1977 and 1995, respectively, was analysed (Table 3). The percent difference between the 
genotype groups did not change significantly fi-om 1977 to 1995, which implied that the 
OPG-A163G polymorphism exerted an influence on peak bone mass. The two 
polymorphisms, called BSP-A1496G (SEQ ID NO. 13) and BSP-G1869A (SEQ ID NO. 14), 
described above also have an impact on peak bone mass. Hence, it was of interest to examine 
whether any co-operation between the OPG polymorphism and either of the BSP 
polymorphisms existed. The combinations OPG-A163G/BSP-A1496G and OPG- 
A163G/BSP-G1869A showed that those polymorphisms certainly act in a co-operative 
fashion, in that the t-test p-values for the null-hypothesis (i.e. no difference between the 
genotype groups) dropped to statistically significant values and the percent difference in mean 
BMC/BMD values for two genotype groups increased (Table 3). 

Please replace the paragraphs beginning at page 38, lines 12 through page 39, line 19 
with the following: 

Thus, it is clear that the OPG-A163G polymorphism, especially in combination with 
the BSP-A1496G (SEQ ID NO. 13) and BSP-G1869A (SEQ ID NO. 14) polymorphisms, is a 
good site for predicting whether an individual is genetically predisposed for high or low 
BMC/BMD. Age, height and weight of the individuals involved in the 18 year study did not 
differ significantly between any of the genotype groups (data not shown). 

To substantiate the initial indication, that the OPG/OCIF polymorphism influences 
peak bone mass, an analysis of the variation of BMC as measured 4 times on the same 
individual from 1977 to 1995 for the different genotypes was performed. In 1977 the 
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average age of the individuals included in the 18 year study was 51.1 years, thus ending at 
69.1 years in 1995. The expressed outcome of a plot of the BMC means for the OPG-A163G 
polymorphism as a function of time would be two parallel curves, each representing BMCs 
measured in individuals with the wild type genotype and BMCs measured in individuals with 
the heterozygous or polymorphic homozygous phenotypes. Figure 1 1 shows that this is 
certainly the case. Moreover, the combination of the OPG-A163G and the BSP-A1496G 
(SEQ ID NO. 13) polymorphisms show that they act in concert on peak bone mass to 
augment the mean HMD difference between genotypes even more than the isolated 
contribution of each polymorphism (Figure 12). In fact, this co-operation is completely 
additive (Table 4), indicating that the two polymorphisms act on bone mass independently of 
one another. The numbers for the BSP-A1496G (SEQ ID NO. 13) and BSP-G1869A (SEQ 
ID NO. 14) polymorphisms in Table 4 are from the results presented in Figures 2 and 3. 
Also, the combination of the OPG-A163G and BSP-G1869A polymorphisms indicates a 
positive co-operation (Figure 13), which is almost additive (Table 4). 

IN THE CLAIMS: 

Please amend the claims as follows: 

4. (Once Amended) A method as claimed in Claim 3, wherein said allelic 
variation of the bone sialoprotein gene promoter is BSP-A1496G (SEQ ID NO. 13) or 
BSP-G1869A (SEQ ID NO 14). 

5. (Once Amended) A method as claimed in Claim 3, wherein said allelic 
variation of the matrix gla protein gene promoter is MGP-C242A (SEQ ID NO. 15). 

6. (Once Amended) A method as claimed in Claim 3, wherein said allelic 
variation of the osteopontin gene promoter is OPN-G520A (SEQ ID NO. 16) or 
OPN-T1825C (SEQ ID NO. 17). 

9. (Once Amended) A method as claimed in Claim 8, wherein said allelic 
variation of the bone sialoprotein gene promoter is BSP-A1496G (SEQ ID NO. 13) or 
BSP-G1869A (SEQ ID NO. 14). 
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IN THE SEQUENCE LISTING: 

Please enter the Sequence Listing submitted herewith, on 3 numbered pages, 
into the instant application in accordance with 32 C.F.R. §§ 1.77(a)(15) and 1.821(c). 
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REMARKS 



In order to satisfy formal requirements, the specification has been amended as shown 
in Appendix A, attached hereto, and the claims have been amended as shown in Appendix B, 
attached hereto. No new matter has been added. 

Pursuant to 37 C.F.R. § 1.821(f), the content of the paper and computer readable 
copies of the Sequence Listing enclosed herewith are the same. No new matter has been 
added. 

Applicant respectfully requests entry of the above-made amendments, the 
accompanying Sequence Listing, and the foregoing remarks into the file of this application. 

It is not believed that fees are required beyond those be provided for in documents 
accompanying this paper. However, if additional fees are required, please charge such fees to 
Pennie & Edmonds LLP Deposit Account No. 16-1 150. 



Respectfully submitted, 



Date: February 5. 2002 




(Reg. No. 33,821) 



For: Paul Zegger 



PENNIE & EDMONDS llp 

1 667 K Street, N.W. 
Washington, D.C. 20006 
(202) 496-4400 



Enclosures 
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APPENDIX A 



Please replace the paragraph beginning at page 17, lines 10 through 23 with the 
following: 

Figure 1 panel A shows the location of the two polymorphisms, called BSP-A1496G 
rSEOIDNO. 13) and BPS-Gl 869A rSEO ID NO. U\ in the bone sialoprotein gene 
promoter. Panel B shows the location of the polymorphism, called MGP-C242A (SEQ ID 
NO. 15\ in the matrix gla protein gene promoter. Panel C shows the location of the 
polymorphisms, called OPN-G520A fSEOIDNO. 16) and OPN-T1825C (SEO ID NO. 17). 
in the osteopontin gene promoter. The wild type sequences encompassing the four 
polymorphic sites for all three said promoters are shown with the nucleotide at the 
polymorphic position in bold and with the substituting nucleotide - also in bold - positioned 
above the polymorphic site. All nucleotide numbering is relative to base pair number 1, 
which is the most 5' nucleotide of each of the promoter sequences as published in the 
GenBank nucleotide database. 

Please replace the paragraphs bridging pages 22 through 24 beginning at page 22, line 
21 through page 24, line 20 with the following: 

DNA analyses. Screening for the BSP-A1496G (^SEOIDNO. 13) and the BSP-Gl 869A 
(SEP ID NO. 14) polymorphisms (basepair numbering according to numbering of BSP 
promoter sequence submitted to GenBank, accession #L24756), the MGP-C242A (SEQ ID 
NO. 15) polymorphism (basepair numbering according to numbering of MGP promoter 
sequence submitted to GenBank, accession #M55270), as well as the OPN-G520A rSEO ID 
NO. 16) and OPN-T1825C rSEO ID NO. 17) polymorphisms (basepair numbering according 
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to numbering of osteopontin promoter sequence submitted to GenBank, accession #D1481 3) 
were performed as follows: 

The polymerase chain reaction (PGR) was used to ampUfy approximately 250 bp long 
DNA fragments of the BSP, MGP, and OPN promoters encompassing the BSP-A1496G 
rSEO ID NO. 13). BSP-G1869A rSEOIDNQ. 14) . MGP-C242A rSEOIDNQ. 15). OPN- 
G520A rSEOIDNO. 16) . and OPN-T1825C rSEOIDNO. 17) polymorphic basepairs. PGR 
techniques are well known in the art and it would be within the ambit of a person of ordinary 
skill in this art to identify primers for amplifying a suitable section of the BSP, MGP and 
OPN genes including the positions 1496bp and 1869pb in the BSP promoter, the position 
242bp in the MGP promoter, and the positions 520bp and 1825bp in the osteopontin 
promoter. PGR techniques are described for example in patents US4683202 or 
EP0200362B1. Two hundred ng of genomic DNA was added to 25 jxl reaction containing Ix 
Taq polymerase buffer with 12.5 mM MgCl2 (Perkin Elmer), 5 nmol of each dNTP, 20 pmol 
of forward and reverse primer, and 1.25 units of AmpliTaq Gold (Perkin Elmer). The 
reaction was heated to 95 °C for 9 minutes followed by 35 cycles of 95 °G for 30 seconds, 
46°G (BSP-A1496G and BSP-G1869A polymorphisms) or 49°C (MGP-C242A 
polymorphism) or 46°G (OPN-G520A polymorphism) or 48°C (OPN-T1825G) for 30 
seconds and 72 °G for 30 seconds - the latter incubation with a 5 second time extension per 
cycle. The reaction was finally incubated 7 minutes at 72°G for completion of the extension 
reaction. Primer sequences for PGR amplification of DNA fragments encompassing the BSP- 
A1496G, BSP-G1869A, MGP-G242A, and OPN-G520A and OPN-T1825G polymorphic 
basepairs were: 

BSP-1496G polymorphism primer set: 
Forward primer: 5' - GAA AAG ATA TAT ATA GAA GGG GAA G - 3' (SEQ ID No. 1) 
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Reverse primer: 5' - TAA TAT CAT TTG ATG TTT CCT CCT G - 3' (SEQ ID No. 2) 

BSP-G1869A polymorphism primer set: 
Forward primer: 5' - TTC TTT CGA CAT AGT GAA AAC ACG T - 3' (SEQ ID No. 3) 
Reverse primer: 5' - CGT GGA TTC TCA CCA GAA AAC - 3' (SEQ ID No. 4) 

MGP-C242A polymorphism primer set: 
Forward primer: 5' - CAG TGA GAA AGC TCA TCA CTT GGT C - 3' (SEQ ID No. 5) 
Reverse primer: 5' - ATT CTC CCA TCC ATC CAT CCA TGC A - 3' (SEQ ID No. 6) 

OPN-G520A polymorphism primer set: 
Forward primer: 5' - CGC TGG AAT TAA GAA AAT TGG TAG A - 3' (SEQ ID No. 7) 
Reverse primer: 5' - GTT GTC AAT TTA GTG GAG GGA GAT C - 3' (SEQ ID No. 8) 

OPN-T1825C polymorphism primer set: 
Forward primer: 5' - GAG TAG TAA AGG ACA GAG GCG AGC T - 3' (SEQ ID No. 9) 
Reverse primer: 5' - CTA GCT TTT TCA TTT ACG GGA TGG G - 3' (SEQ ID No. 10) 

Please replace the paragraph begimiing at page 27, lines 15 through 27 with the 
following: 

Results 

Five previously unknown polymorphisms were identified by sequencing specific 
promoter regions from the human BSP gene promoter, the hiiman MGP gene promoter, and 
the human OPN gene promoter following a PCR amplification of 40 DNA samples from 
healthy women. The BSP-Al 496G (SEP ID NO. 131 BSP-Gl 869A (SEP ID NO. 141 
MGP-C242A rSEOIDNO. 15\ OPN-G520A (SEP ID NO. 16) and OPN-Tl 825C rSEO ID 
NO. 17) polymorphisms were coded as Xx, Yy, Zz, Bb, and Ss, respectively, where the 
uppercase letter signifies presence of the wild type base pair at the given poljonorphic 



position and the lowercase letter signifies presence of the base pair different from the wild 
type base pair at the given polymorphic position. 

Please replace the paragraph bridging pages 27 and 28 beginning at page 27, lines 30 
through page 28, line 9 with the following: 

Table 1 shows the genotype distribution of DNA samples from the 18 years study for 
all identified polymorphic sites: BSP-A1496G rSEO ID NO. 13) . BSP-Gl 869A rSEO ID 
NO. 14). MGP-C242A (SEP ID NO. 15) . OPN-G520A fSEO ID NO. 16) and OPN-T1825C 
(SEP ID NO. 17) . The left panel shows the actual number of samples categorized into three 
genotypes for the 3 identified polymorphic sites. The right panel displays the same analysis 
as the left except that the numbers represent the percent of total samples analyzed for each 
polymorphic site. 

wt = XX, YY, ZZ, BB or SS 

hz = Xx, Yy, Zz, Bb or Ss 

pm = XX, yy, zz, bb or ss 

Please replace the paragraph beginning at page 29, lines 3 through 20 with the 
following: 

The genotype distributions for the 5 polymorphisms are shown in table 1 . For the 
BSP-A1496G (SEP ID NP. 13) polymorphism the homozygous wild type genotype was the 
most abundant, followed by the heterozygous and homozygous polymorphic genotypes, with 
the homozygous polymorphic groups being quite small. In the case of BSP-Gl 869A (SEP 
ID NO. 14) polymorphism, the wild type genotype, as defined by the BSP gene promoter 
sequence from GenBank, was rare, the heterozygous genotype was 10 times more frequent 
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and the homozygous polymorphic genotype was twice as frequent as the heterozygous. For 
the MGP-C242 (SEP ID NO. 15) polymorphism the heterozygous genotype was the most 
abundant followed by the homozygous wild type and homozygous polymorphic genotypes. 
In the case of the OPN-G520A rSEOIDNO. 16) polymorphism, the homozygous 
polymorphic genotype was the most abundant followed by the heterozygous and the wild type 
homozygous genotypes. The genotype distribution of the OPN-T1825C (SEP ID NO. 17) 
polymorphism was, generally, the same as for the MGP-C242A fSEOIDNO. 15) 
polymorphism. 

Please replace the paragraphs beginning at page 32, lines 1 through page 34, line 16 
with the following: 

From the table it is clear that the BSP-A1496G CSEOIDNO. 13) and BSP-Gl 869A 
(SEP ID NO. 14) polymorphic sites, especially when combined, are good sites for predicting 
whether an individual is genetically predisposed for high or low BMC/BMD. The PPN- 
Tl 825C (SEP ID NO. 17) polymorphism only has a marginal influence on BMC/BMD on its 
own. However, when the PPN-T1825C (SEPIDNP. 17) polymorphism is combined with 
the BSP-Gl 869A (SEP IDNP. 14) polymorphism the percent separation of genotypes is 
better than either polymorphism alone (Table 2). Pn the other hand, the MGP-C242A (SEP 
IDNP. 15) and PPN-G520A (SEP IDNP. 16) polymorphisms are, at first glance, not 
suitable sites for such a prediction. None of the identified polymorphisms appeared to have a 
statistically significant impact on the change in bone mass over time (data not shown). Age, 
height and weight of the individuals involved in the 18 years study did not differ significantly 
between any of the genotype groups (data not shown). 
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These observations strongly indicate that the BSP polymorphisms influence peak bone 
mass rather than the rate of bone loss. To substantiate this an analysis of the variation of 
BMD as measured 4 times on the same individual from 1977 to 1995 for the different 
genotypes was performed. In 1977 the average age of the individuals included in the 18 years 
study v^as 51.1 years, thus ending at 69.1 years in 1995. The expected outcome of a plot of 
the means of BMD for one of the two BSP polymorphisms as a function of time would be 
two parallel curves, each representing BMDs measured in individuals vdth the wild type 
genotype and BMDs measured in individuals with the polymorphic phenotype. Figures 2 and 
3 show that this is, indeed, the case for the BSP-A1496G (SEP ID NO. 13) and BSP-Gl 869A 
(SEP ID NO. 14) polymorphic sites. Moreover, the two BSP promoter polymorphism act in 
concert on peak bone mass to augment the mean BMD difference between genotypes even 
more than the isolated contributions of each polymorphism (Figure 4). 

The role - if any - of the MGP-C242A (SEPIDNP. 15) and PPN-G520A (SEP ID 
NP. 16) polymorphisms in the MGP and PPN promoters on bone turnover was less clear 
from the first analyses compiled in table 2. However, when BMC values grouped according 
to genotype were plotted as a function of time a set of curves appeared suggesting that both 
the MGP-C242A rSEPIDNP. 15) (Figure 5) and OPN-G520A (SEP ID NO. 16) (Figure 6) 
polymorphic sites are determinants of rate of bone loss. It is especially noteworthy that the 
ZZ and Zz+zz curves as well as the BB+Bb and bb curves separate between 1979 and 1989, 
corresponding to an average age of 53.1 years and 63.1 years, indicative of a genetic 
phenomenon associated with the menopause. Like the BSP polymorphisms, the combined 
action of the MGP-C242A (SEPIDNP. 15) and PPN-G520A (SEP IDNP. 16) 
polymorphisms also leads to a bigger difference between genotypes than either would create 
alone (Figure 7). 
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According to the results compiled in table 2, the impact of the OPN-T1825C fSEQ ID 
NO. 11) polymorphism on BMC/BMD was only visible after it was combined with the BSP- 
G1869A (SEP ID NO. 14) polymorphism. From a graph of BMC values grouped according 
to genotype and plotted as a function of time it is difficult to tell whether this polymorphism 
has an impact on rate of bone loss or peak bone mass, due to the proximity of the curves 
(Figure 8). However, combining this polymorphism with the BSP-Gl 869 A fSEO ID NO. 14) 
polymorphism gave rise to a set of time course curves clearly showing that these two 
polymorphisms cooperate in an additive fashion, and, hence, that the OPN-T1825C rSEO ID 
NO. 17) polymorphism may influence peak bone mass, as the BSP-Gl 869A (SEP ID NO. 
14) polymorphism, rather than the rate of bone loss (Figure 9). 

Finally, the association between genotype and urinary osteocalcin (N-MID®), 
Psteometer Biotech A/S) as well as urinary collagen type 1 C-terminal crosslinks 
(CrossLaps® Osteometer Biotech A/S) were examined. As expected there was no significant 
difference between the mean value of either of the biochemical bone turnover markers for the 
BSP-A1496G (SEPIDNP. 13) and BSP-Gl 869A (SEP ID NO. 14) polymorphic sites (data 
not shown). Also, no significant difference was observed between the mean value of the 
biochemical bone turnover markers for the MGP-C242A fSEOIDNO. 15) and OPN-G520 A 
(SEP ID NO. 16) polymorphic sites. This is likely due to the equally paced bone loss of the 
ZZ and Zz+zz genotype groups as well as the BB+Bb and bb genotypes groups at the time of 
N-MID® and CrossLaps® measurement (1995) according to Figures 5 and 6. 

Please replace the paragraph beginning at page 37, lines 1 through 19 with the 
following: 
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The impact of the identified polymorphic site on bone mass as represented by bone 
mineral content (BMC) and bone mineral density (BMD) measurements at the distal arm in 
1977 and 1995, respectively, was analysed (Table 3). The percent difference between the 
genotype groups did not change significantly from 1977 to 1995, which implied that the 
OPG-A163G polymorphism exerted an influence on peak bone mass. The two 
polymorphisms, called BSP-A1496G (SEP ID NO. 13^ and BSP-G1869A rSEOIDNO. U\ 
described above also have an impact on peak bone mass. Hence, it was of interest to examine 
whether any co-operation between the OPG polymorphism and either of the BSP 
polymorphisms existed. The combmations OPG-A163G/BSP-A1496G and OPG- 
A163G/BSP-G1869A showed that those [polymor-phisms] polymorphisms certainly act in a 
co-operative fashion, in that the t-test p-values for the null-hypothesis (i.e. no difference 
between the genotype groups) dropped to statistically significant values and the percent 
difference in mean BMC/BMD values for two genotype groups increased (Table 3). 

Please replace the paragraphs beginning at page 38, lines 12 through page 39, line 19 
with the following: 

Thus, it is clear that the OPG-Al 63G polymorphism, especially in combination with 
the BSP-A1496G (SEP ID NO. 13^ and BSP-G1869A CSEO ID NO. W) polymorphisms, is a 
good site for predicting whether an individual is genetically predisposed for high or low 
BMC/BMD. Age, height and weight of the individuals involved in the 18 year study did not 
differ significantly between any of the genotype groups (data not shown). 

To substantiate the initial indication, that the OPG/OCIF polymorphism influences 
peak bone mass, an analysis of the variation of BMC as measured 4 times on the same 
individual from 1 977 to 1 995 for the different genotypes was performed. In 1 977 the 
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average age of the individuals included in the 18 year study was 51.1 years, thus ending at 
69.1 years in 1995. The expressed outcome of a plot of the BMC means for the OPG-A163G 
polymorphism as a function of time would be two parallel curves, each representing BMCs 
measured in individuals with the wild type genotype and BMCs measured in individuals with 
the heterozygous or polymorphic homozygous phenotypes. Figure 1 1 shows that this is 
certainly the case. Moreover, the combination of the OPG-A163G and the BSP-A1496G 
rSEO ID NO. 13) polymorphisms show that they act in concert on peak bone mass to 
augment the mean BMD difference between genotypes even more than the isolated 
contribution of each polymorphism (Figure 12). In fact, this co-operation is completely 
additive (Table 4), indicating that the two polymorphisms act on bone mass independently of 
one another. The numbers for the BSP-A1496G rSEO ID NO. 13) and BSP-Gl 869A (SEP 
ID NO. 14) polymorphisms in Table 4 are from the results presented in Figures 2 and 3. 
Also, the combination of the OPG-A163G and BSP-Gl 869A polymorphisms indicates a 
positive co-operation (Figure 13), which is almost additive (Table 4). 
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APPENDIX B 



Pending claims 
Please amend the claims as follows: 

4. (Once Amended) A method as claimed in Claim 3, wherein said allelic 
variation of the bone sialoprotein gene promoter is BSP-A1496G (SEP ID NO. 13) or 
BSP-Gl 869A rSEO ID NO 141 

5. (Once Amended) A method as claimed in Claim 3, wherein said allelic 
variation of the matrix gla protein gene promoter is MGP-C242A (SEOIDNO. \5\ 

6. (Once Amended) A method as claimed in Claim 3, wherein said allelic 
variation of the osteopontin gene promoter is OPN-G520A fSEO ID NO. 16) or 
OPN-T1825C rSEOIDNO. 171 

9. (Once Amended) A method as claimed in Claim 8, wherein said allelic 
variation of the bone sialoprotein gene promoter is BSP-A1496G fSEO ID NO. 13) or 
BSP-Gl 869A (SEOIDNO. 141 
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^- ' 09/889491 

JCt8R8C'dPCT/PT0 I 8 JUL 2001 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

Application of: P. HUSK 

National Stage of PCT/EPOO/003 1 9 Group Art Unit: Unassigned 

Filed: July 18, 2001 Examiner: Unassigned 

For: GENETIC PREDISPOSITION Attorney Docket No. : 8969-029 

PRELIMINARY AMENDMENT 

Box PATENT APPLICATION 

Assistant Commissioner for Patents 
Washington, D.C. 20231 

Sir: 

Please enter the following amendments and remarks into the file of the above- 
identified application prior to the examination thereof. 

IN THE SPECIFICATION 

Marked up versions of all paragraphs showing insertions and deletions are 
included in Appendix A. 

Please replace the paragraph starting at page 1, line 5 with the following text: 
-TECHNICAL FIELD 

This invention relates to a method for assessing predisposition to various conditions 
based upon polymorphisms in a bone sialoprotein gene, a matrix gla protein gene, an 
osteopontin gene and/or an osteoprotegerin (OPG)/osteoclastogenesis inhibitory factor 
(OCIF) gene. More specifically, the invention relates to a method of assessing an 
individual's predisposition to various pathological calcification conditions including osteo- 
porosis and atherosclerosis by screening for these polymorphisms. The method of the present 
invention is especially useful in determining allehc variations in the himian bone sialoprotein 
gene, the human matrix gla protein gene, the himian osteopontin gene and/or the 
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osteoprotegerin (OPG) /osteoclastogenesis inhibitory factor (OCIF) gene thus predicting 
predisposition to high or low bone mineral density (BMD). The invention also relates to 
bone sialoprotein (BSP) genes, matrix gla protein (MGP) genes, ostepeopontin (OPN) genes 
and OPG/OCIF genes containing the polymorphisms and to probes and primers therefor.-- 

Please replace the paragraph starting at page 5, Une 16 with the following text: 
--SUMMARY OF THE INVENTION 

The present invention now provides a method of assessing an individual's 
predisposition to a selected calcification condition status, which method comprises 
determining the genotype of the promoter of the bone sialoprotein gene, the promoter of the 
matrix gla protein gene, the promoter of the osteopontin gene, or the promoter of the 
OPG/OCIF gene or all four or any combination of two or more out of the four promoters.-- 

Please replace the paragraph starting at page 17, line 6 with the following text: 
-BRIEF DESCRIPTION OF THE DRAWINGS 

In order that the nature of the present invention be more clearly understood, there now 
follows an example in which reference is made to the Figures shown in the accompanying 
drawings in which: — 

Please replace the paragraph starting at page 22, line 13 with the following text: 
-DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Example 1 (18 year study) 

Methods 

Subjects. One hundred thirty three women followed up for 18 years (1977-1995) with respect 
to BMD, biochemical markers, height, and weight were used in the study. A detailed 
description of the cohort has been previously published (Jorgensen et al., 1996).— 

IN THE CLAIMS 

Please rewrite claims 3, 14, 23, and 24 to recite.: 

A list of changes to the claims showing insertions and deletions is provided herewith 
as Appendix B. 
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4. - 

3. (Amended once) A method as claimed in Claim 1 , wherein it is determined 
whether the individual is homozygous or heterozygous for an allelic variation of the promoter 
of the bone sialoprotein gene, the promoter of the matrix gla protein gene, the promoter of the 
osteopontin gene or the promoter of the OPG/QCIF gene, or all four or a combination of two 
or more out of the four promoters. 

14. ( Amended once) A method as claimed in Claim 1 , comprising amphfying a 
relevant portion of the DNA of a said gene promoter of said individual. 

23. (Amended once) A method of osteoporosis therapy comprising determining a 
predisposition as claimed in Claim 1, and administering a medicament to the individual to 
prevent or treat osteoporosis or to delay the onset of osteoporosis if the individual is 
predisposed to low peak bone mass or to a high rate of loss of bone mass. 

24. (Amended once) A method of atherosclerosis therapy comprising determining 
a predisposition as claimed in Claim 1, and administering a medicament to the individual to 
prevent or treat atherosclerosis or to delay the onset of atherosclerosis if the individual is 
predisposed to pathological arterial calcification. 
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REMARKS 

Claims 1-24, as amended, are pending in this application for the Examiner's 
review and consideration. Applicants have amended the specification and claims to conform 
with U.S. patent practice and to more clearly recite the invention. Claims 3, 14, 23, and 24 
were amended so that they are not multiple dependent claims. As no new matter has been 
added herein, these changes should be entered. 



Respectfully submitted. 

Date yJL^ 2do\ f ^ SWi^^T". L-^^-^^ ^ TH T 

^ Wanton T. Lawrence, IE (Reg. No. 25,736) 



PENNIE & EDMONDS llp 

1667 K Street, N.W. 
Washington, DC 20006 

(202) 496-4460 
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Appendix A 

Changes to the Specification 
The paragraph starting at page 1, hne 5 is amended as follows: 
- TECHNICAL FIELD 

This invention relates to a method for assessing predisposition to various conditions 
based upon polymorphisms in a bone sialoprotein gene, a matrix gla protein gene, an 
osteopontin gene and/or an osteoprotegerin (OPG)/osteoclastogenesis inhibitory factor 
(OCIF) gene. More specifically, the invention relates to a method of assessing an 
individual's predisposition to various pathological calcification conditions including osteo- 
porosis and atherosclerosis by screening for these polymorphisms. The method of the present 
invention is especially useful in determining allelic variations in the human bone sialoprotein 
gene, the human matrix gla protein gene, the human osteopontin gene and/or the 
osteoprotegerin (OPG) /osteoclastogenesis inhibitory factor (OCIF) gene thus predicting 
predisposition to high or low bone mineral density (BMD). The invention also relates to 
bone sialoprotein (BSP) genes, matrix gla protein (MGP) genes, ostepeopontin (OPN) genes 
and OPG/OCIF genes containing the polymorphisms and to probes and primers therefor. -- 

The paragraph starting at page 5, line 16 is amended as follows: 
- SUMMARY OF THE INVENTION 

The present invention now provides a method of assessing an individual's 
predisposition to a selected calcification condition status, which method comprises 
determining the genotype of the promoter of the bone sialoprotein gene, the promoter of the 
matrix gla protein gene, the promoter of the osteopontin gene, or the promoter of the 
OPG/OCIF gene or all four or any combination of two or more out of the four promoters. - 

The paragraph starting at page 17, line 6 is amended as follows: 
- BRIEF DESCRIPTION OF THE DRAWINGS 

In order that the nature of the present invention be more clearly imderstood, there now 
follows an example in which reference is made to the Figures shown in the accompanying 
drawings in which: — 
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The paragraph starting at page 22, line 13 is amended as follows: 
-DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Example 1 (18 year study) 

Methods 

Subjects. One hundred thirty three women followed up for 18 years (1977-1995) with respect 
to BMD, biochemical markers, height, and weight were used in the study. A detailed 
description of the cohort has been previously published (Jorgensen et al., 1996).— 
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Appendix B 

Changes to the Claims 



3. (Amended once) A method as claimed in Claim 1 [or Claim 2], wherein it is 
determined whether the individual is homozygous or heterozygous for an allelic variation of 
the promoter of the bone sialoprotein gene, the promoter of the matrix gla protein gene, the 
promoter of the osteopontin gene or the promoter of the OPG/QCIF gene, or all fom or a 
combination of two or more out of the four promoters. 

14. ( Amended once) A method as claimed in Claim 1 [any preceding claim], 
comprising amplifying a relevant portion of the DNA of a said gene promoter of said 
individual. 

23. (Amended once) A method of osteoporosis therapy comprising determining a 
predisposition as claimed in [any one of Claims 1 to 22] Claim 1 , and administering a 
medicament to the individual to prevent or treat osteoporosis or to delay the onset of 
osteoporosis if the individual is predisposed to low peak bone mass or to a high rate of loss of 
bone mass. 

24. (Amended once) A method of atherosclerosis therapy comprising determining 
a predisposition as claimed in [any one of Claims 1 to 22] Claim 1, and administering a 
medicament to the individual to prevent or treat atherosclerosis or to delay the onset of 
atherosclerosis if the individual is predisposed to pathological arterial calcification. 
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09/886491 
WWPCTR9C*«J03FEB2002 

SEQUENCE LISTING 

<110> Philip KUSK 

<120> Genetic Predisposition 

<130> 8969-029-999 



<140> 09/889,491 
<141> 2001-07-18 



<150> 9901037.3 
<151> 1999-01-18 



<150> 9912585.8 
<151> 1999-05-28 



<150> PCT/EPOO/00319 
<151> 2000-01-17 



<160> 17 



<170> FastSEQ for Windows Version 4.0 

<210> 1 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 1 

gaaaagatat atatagaagc ccaag 

<210> 2 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 2 

taatatcatt tgatgtttcc tcctg 

<210> 3 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 3 

ttctttcgac atagtgaaaa cacgt 

<210> 4 

<211> 21 

<212> DNA 

<213> Homo sapiens 

<400> 4 

cgtggattct caccagaaaa c 



<210> 5 

<211> 25 

<212> DNA 

<213> Homo sapiens 



_1. 



<400> 5 

cagtgagaaa gctcatcact tggtc 



<210> 6 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 6 

attctcccat ccatccatcc atgca 

<210> 7 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 7 

cgctggaatt aagaaaattg gtaga 

<210> 8 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 8 

gttgtcaatt tagtggaggg agate 

<210> 9 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 9 

gagtagtaaa ggacagaggc gagct 

<210> 10 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 10 

ctagcttttt catttacggg atggg 

<210> 11 

<211> 25 

<212> DNA 

<213> Homo sapiens 

<400> 11 

agtctaactt ctagaccagg caatt 

<210> 12 

<211> 19 

<212> DNA 

<213> Homo sapiens 

<400> 12 

agttagagca gagaatctg 



<210> 13 

<211> 30 

<212> DNA 

<213> Homo sapiens 

<400> 13 

atatagaagc ccaagaaaaa tcagctgacc 

<210> 14 

<211> 31 

<212> DNA 

<213> Homo sapiens 

<400> 14 

atagtgaaaa cttgtgtaat tatgaaattt 

<210> 15 

<211> 26 

<212> DNA 

<213> Homo sapiens 

<400> 15 

tggctggctg gctggctgga tggatg 

<210> 16 

<211> 30 

<212> DNA 

<213> Homo sapiens 

<400> 16 

gtttttagaa ttttcggact tccctccact 

<210> 17 

<211> 30 

<212> DNA 

<213> Homo sapiens 



<400> 17 

gacagaggca agttttctga agtccttgca 
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GENETIC PREDISPOSITION 



This invention relates to a method for assessing 
predisposition to various conditions based upon 
polymorphisms in a bone sialoprotein gene, a matrix gla 
protein gene, an osteopontin gene and/or an osteoprotegerin 
(OPG) /osteoclastogenesis inhibitory factor (OCIF) gene. 
More specifically, the invention relates to a method of 
assessing .an individual's predisposition to various 
pathological calcification conditions including osteo- 
porosis and atherosclerosis by screening for these 
polymorphisms. The method of the present invention is 
especially useful in determining allelic variations in the 
human bone sialoprotein gene, the human matrix gla protein 
gene, the human osteopontin gene and/or the osteoprotegerin 
{OPG) /osteoclastogenesis inhibitory factor (OCIF) gene thus 
predicting predisposition to high or low bone mineral 
density (BMD) . The invention also relates to bone 
sialoprotein (BSP) genes, matrix gla protein (MOP) genes, 
ostepeopontin (OPN) genes and OPG/OCIF genes containing the 
polymorphisms and to probes and primers therefor. 

Osteoporosis is today one of the most common diseases 
in individuals over 60 years of age. In America alone it 
affects an estimated 25 million people with a 5:1 ratio of 
women to men. This corresponds to approximately 25-3 0% of 
people over 60 years of age. In Europe the percentage of 
people affected by this disease is approximately the same. 

At present there is no cure for osteoporosis. However, 
hormone replacement therapy as well as treatment with 
bisphosphonates can halt or slow accelerated bone loss. 
Hence, the sooner such a bone loss can be diagnosed the 
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better the impact of treatment. It would, accordingly, be 
of particular advantage to be able to identify individuals 
predisposed to osteoporosis as early as possible. Since 
there is a strong genetic component associated with the 
5 development of osteoporosis (see below) the identification 
of genes having an impact on peak bone mass and/or the rate 
of bone loss would be of enormous help because individuals 
with genotypes predisposing to osteoporosis could be 
identified early in their life, leaving plenty of time for 

10 preventive .measures to be instituted. 

From this point of view it is relevant to identify 
individuals or groups who are predisposed either to have a 
relatively low peak bone mass or bone mineral density or 
who are predisposed to a relatively fast rate of loss of 

15 bone mass or bone mineral density. 

The extent of genetic contribution to bone mineral 
density and, hence, possibly to osteoporosis became evident 
more than 20 years ago from twin studies (Smith et al . , 
1973) . Larger BMD and bone width variances were found in 

20 dizygotic than in monozygotic twins, indicating a 
considerable genetic influence on the regulation of bone 
mass (Smith et al . , 1973). Another twin study reported a 
statistically significant effect of genetic factors on the 
rate of BMD change (Kelly et al . , 1993). Other studies have 

25 demonstrated a strong genetic impact on the acquisition of 
peak bone mass (Gueguen et al . , 1995,- Lutz & Tesar, 1990). 
Hence, the genetic component seemingly has an impact on 
both peak bone mass and the rate of bone loss . Genetic 
segregation analysis has strongly suggested that bone mass 

30 is controlled by several genes, each with modest effects 
(Gueguen et al . , 1995). 
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Very little is known about the molecular mechanisms 
and genetics, leading to osteoporosis, even though there 
has been an intense search for genes influencing bone 
formation and bone resorption. Within the last five years 
polymorphisms identified in genes mainly encoding 
regulatory proteins have been associated with BMD: 

1) A vitamin D receptor gene polymorphism identified by 
Morrison and co-workers (Morrison et al . , 1994; 
WO94/03633) has been subject of numerous publications. 
Only about half found a significant but weak association 
and this association may have a stronger link to peak 
bone mass than rate of bone loss as suggested by Riggs 
(Riggs, 1997) . 

2) An estrogen receptor gene polymorphism identified by 
Kobayashi and co-workers has been associated with low BMD 
(Kobayashi et al . , 1996, US5834200). 

3) A polymorphism in the collagen type I a 1 gene 
identified by Grant and co-workers (Grant et al . , 1996; 
W09732041) has been associated with rate of bone loss 
rather than peak bone mass (Uitterlinden et al . , 1998) . 

4) Two polymorphisms in an AT rich minisatellite repeat 
in the 3' flank of the interleukin- 6 gene - an 
interleukin which has a stimulatory effect on cells of 
the osteoclast lineage - have been found to be associated 
with low BMD (Murray et al . , 1997; W09743446) . These 
polymorphisms seem to influence peak bone mass. 

5) Two tandem repeat polymorphisms recently found in the 
interleukin- 1 receptor antagonist gene have been 
associated with rate of bone loss (Keen et al . , 1998; 
WO9844150) . The reason for focusing on this gene was that 
interleukin- 1 is a potent stimulator of osteoclastic bone 
resorption. Interleukin- 1 is stimulated and interleukin- 1 
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receptor antagonist is inhibited in estrogen deficient 
states, where the associated bone loss can be blocked by 
treatment with interleukin- 1 receptor antagonist, which 
then seemed to be a good candidate gene for the 
5 regulation of postmenopausal bone loss (Keen et al . , 

1998) . 

6) A polymorphism in the osteocalcin gene has recently 
been reported {Dohi et al . , 1998). This polymorphism was 
unable to discriminate BMDs measured in persons with the 

10 aberrant .allele from BMDs measured in persons with the 

wild type genotype in a statistically significant manner. 

7) Polymorphisms linked to low HMD or osteoporosis have 
been described in the TGF-Sl gene, whose protein product 
is abundant in bone and an important regulator of bone 

15 resorption and formation (Langdahl et al . , 1997; Yamada 
et al., 1998; WO97/28280) . 

8) A polymorphism in the apolipoprotein E gene has been 
reported to be correlated to low BMD in postmenopausal 
Japanese women {Shiraki et al., 1997). The reason for 

20 focusing on this gene was a report stating that the level 
of vitamin K which activates osteocalcin through y- 
glutarayl carboxylase was related to the apolipoprotein E 
phenotype (Saupe et al , , 1993). 

9) WO9705275 discloses use of analysing for an allelic 
25 variant in the retinoic acid receptor gene for the 

prediction of bone density. 

The vitamin D receptor gene and the collagen type I 
1 gene polymorphisms only seem to be associated with BMD in 
a fraction of the cohorts examined (the other gene 
30 polymorphisms have only been tested on small, national 
cohorts and most of them are rather peripheral to bone 
metabolism) . Since the genetic influence on the development 
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of osteoporosis is caused by the inadequate action of 
multiple genes, this comes as no surprise. Obviously, more 
genetic polymorphisms with an impact on bone formation/ 
resorption need to be identified to get a better genetic 
5 prediction power valid within a wider geographical area. In 
this respect, we believe it is important to focus on the 
promoter regions of genes encoding proteins present in bone 
for the following reason: Mutations predisposing for 
osteoporosis must be rather subtle in nature, since they 

10 generally resist physical detection for more than 50 years. 
Mutations affecting the expression of a gene product in 
stead of its function would be expected to display such 
subtle effects. Gene expression is controlled by promoters 
residing upstream of the coding region of a gene. Mutations 

15 in these regions could cause an altered gene expression. 

The present invention now provides a method of 
assessing an individual's predisposition to a selected 
calcification condition status, which method comprises 
determining the genotype of the promoter of the bone 

20 sialoprotein gene, the promoter of the matrix gla protein 
gene, the promoter of the osteopontin gene, of the promoter 
of the OPG/OCIF gene or all four or any combination of two 
or more out of the four promoters. 

The calcification condition status for which a 

25 predisposition is assessed according to the invention may 
be having a high or a low peak bone mass {as a future, 
present or as a past state) or having a high rate or a low 
rate of bone loss (as a future, present or past state) . 
Thus, the invention may be used to assess a predispost ion 

30 to osteoporosis. 

In a second aspect the present invention provides a 
method of assessing predisposition of an individual to any 
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condition associated with allelic variation of a said 
promoter or any such combination thereof. 

The method of the invention typically comprises 
determining whether an individual is homozygous or heter- 
5 ozygous for a bone sialoprotein promoter (BSP) , a matrix 
gla protein promoter (MGP) , an osteopontin promoter (OPN) , 
or an OPG/OCIF promoter or all four promoters or 
combinations of two or more out of these four promoters and 
particular polymorphisms thereof. The method is 

10 conveniently used to screen for an individual at risk of a 
condition or disease correlated with aberrant production of 
bone sialoprotein, matrix gla protein, osteopontin, or 
OPG/OCIF or all four or two or more out of the four such as 
osteoporosis or atherosclerosis. 

15 A DNA sequence of the human bone sialoprotein promoter 

is known and has been published by Kim, R.H. et al . in 
Matrix Biol. 14: 31-40 (1994). The sequence submitted to 
GenBank by this group with accession # L24756 is referred 
to hereafter as the wild type sequence or the published 

20 sequence. A DNA sequence of the human matrix gla protein 
promoter is known and has been published by C'ancela, L. et 
al. in J. Biol. Chem. 265 (25): 15040-15048 (1990). The 
sequence submitted to GenBank by this group with accession 

# M55270 is referred to hereafter as the wild type sequence 
25 or the published sequence. A DNA sequence of the human 

osteopontin promoter is known and has been published by 
Hijiya et al . in Biochem J., 303: 255-262 (1994). The 
sequence submitted to GenBank by this group with accession 

# D14813 is referred to hereafter as the wild type sequence 
30 or the published sequence. A DNA sequence of the human 

OPG/OCIF promoter is known and has been published by 
Morinaga et al . , Eur. J. Biochem. 254 ( 3 ) : 658 - 69 1 (1998). 
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The sequence submitted to GenBank by this group with 
accession #AB008821 is referred to hereafter as the wild 
type sequence or the published sequence. This terminology 
is not intended to imply that any of these published 
sequences is more prevalent in the population than 
variations thereof or that each or any of them is 
associated with the minimum risk of pathology. The method 
of the invention includes determining whether the 
individual being tested has a bone sialoprotein promoter, a 
matrix gla protein promoter, an osteopontin. promoter, or an 
OPG/OCIF promoter or all four or combinations of two or 
more out of these four promoters which are identical with 
the published sequences (or are identical at selected 
regions of said sequences) or whether that individual has a 
bone sialoprotein promoter, a matrix gla protein promoter, 
an osteopontin promoter, or an OPG/OCIF promoter or all 
four or combinations of two or more out of these four which 
differ from the published sequences (or which differ at 
said selected locations), i.e. are polymorphisms of the 
published sequences, whether homozygous or heterozygous. 

The invention includes a method as described above in 
which one determines the sequence at location 14 96 bp of 
the BSP, in particular whether the sequence at this 
location is A (published) or G, and/or at the location 1869 
bp, and in particular whether the sequence at said location 
is G (published) or A. The locations identified above are 
numbered from the start of the published sequence. In an 
alternative numbering scheme, these locations are -68 3 bp 
and -310 bp from the start of the transcribed sequence of 
the gene. Hereafter, these specific allelic variations are 
indicated using the terminology BSP-A1496G and BSP-G1869A. 
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Similarly, the invention includes such a method in 
which said allelic variation is at MGP location 242 bp 
(numbered from the start of the published sequence or -3157 
bp from the start of the transcribed sequence, and in 
particular whether the sequence at said location reads C 
(published) or A, referred to hereafter as MGP-C242A. 

Similarly, the invention includes such a method m 
which said allelic variation is at OPN location 520 bp 
(numbered from the start of the published sequence or -1748 
bp from the start of the transcribed sequence) and in 
particular whether the sequence at said location reads G 
(published) or A, referred to hereafter as OPN-G520A. 

Similarly, the invention includes such a method in 
which said allelic variation is at OPN location 1825 bp 
(numbered from the start of the published sequence or -443 
bp from the start of the transcribed sequence) and in 
particular whether the sequence at said location reads T 
(published) or C, referred to hereafter as OPN-T1825C. 

The invention includes a method as described above in 
which one determines the sequence at location 163 bp of the 
OPG/OCIF promoter, in particular whether the sequence at 
this location is A (published) or G. The location 
specified above is numbered from the start of the published 
sequence. In an alternative numbering scheme, this 

location is -943 bp from the start of the transcribed 
sequence of the gene. Hereafter, this specific allelic 
variation is indicated using the terminology OPG-A163G. 

We have found that if at least one copy of the BSP 
gene of the individual has G rather than A at position 14 96 
or A rather than G at position 1869, this is associated 
with higher peak bone mass. 
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Similarly, we have found that if at least one copy of 
the MGP gene of the individual has A rather than C at 
position 242, this is associated with a higher rate of loss 
of bone mass . 

5 We have further found that having at least one copy of 

A rather than G at OPN position 52 0 bp is associated with a 
higher rate of loss of bone. 

We have further found that having at least one copy of 
T rather than C at OPN position 1825 bp is associated with 
10 a lower bone mass. 

We have found that if at least one copy of the 
OPG/OCIF gene of an individual has A rather than G at 
position 163 bp, this is associated with higher bone mass. 

The relevant determinations of gene promoter sequences 
15 can be carried out by generally known methods, which 
generally involve amplifying a relevant portion of the DNA 
of a said gene promoter of said individual . The sequence 
of said amplified portion may be determined by 
hybridisation assay or by restriction fragment length 
20 analysis. 

In particular, amplification may be conducted using a 
promoter chosen such that if a selected one of the 
published sequences or the variation of the published 
sequence is present, amplification will produce a new site 

25 at which the amplicon will be cut by a restriction enzyme. 
A different number of restriction fragments will thus be 
produced by enzyme treatment of the amplicon. The 
invention includes an oligonucleotide primer for use in 
amplification of a relevant portion of a said gene 

30 promoter. In particular, the invention includes such a 
promoter selected so as to produce a different restriction 
pattern depending on the presence or absence of a selected 
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variation. Suitable promoters according to the invention 
are described in the examples hereafter. 

The invention includes a method of osteoporosis 
therapy comprising determining a predisposition as 
described above, and administering a medicament to the 
individual to prevent or treat osteoporosis or to delay the 
onset of osteoporosis if the individual is predisposed to 
low peak bone mass or to a high rate of loss of bone mass. 

Bone sialoprotein is a bone tissue specific 3 3-34 kDa 
nascent protein which is extensively modified post- 
translationally by glycosylation, phosphorylation and 
sulfation leading to a final MW of 57 kDa (Oldberg et al . , 
1988; Ecarot-Charrier et al . , 1989; Fisher et al . , 1990; 
Zhang et • al . , 1990). Together with osteopontin, BSP is the 
most abundant non- collagenous protein in the bone matrix 
(Nagata et al . , 1991). It contains an RGD motif that 
mediates cell attachment via integrin class of cell 

surface receptor found on osteoclasts (Oldberg et al.l988; 
Flores et al . , 1992; Ross et al . , 1993). It also has 
several segments of poly-glutamic acid that create a potent 
hydroxyapatite nucleating domain (Hunter Sc Goldberg, 1993) . 
Temporal studies on the localization of endogenous BSP in 
rats (Chen et al . , 1991; Chen et al . , 1992) and in 
transgenic mice {Chen et al . , 1996) have shown, that 
highest expression of BSP occurs in neonatal bones, with 
expression decreasing profoundly with subsequent growth and 
development. BSP has been localized ahead of the 
mineralization front using immune -histochemistry tech- 
niques, suggesting that it is necessary for the initiation 
of bone mineralization (Roach, 1994) . Accordingly, mRNA has 
only been detected in differentiated osteoblasts, 
odontoblasts, and cementoblasts at sites of de novo 
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mineralized tissue formation (Chen et al . , 1991; Chen et 
al . , 1992) . In transgenic mice, with the transgene being an 
appr. 2.7 kb promoter region from rat BSP fused to the 
lucif erase gene, the same expression pattern is observed, 
5 indicating that the appr. 2.7 kb region of the rat promoter 
is sufficient to mediate the bone tissue specific 
transcription (Chen et al., 1996) . 

Matrix gla protein is a small 79 amino acid residues 
protein with molecular weight of appr. 14 kDa which 

10 contains fi.ve y-carboxyglutamic acid (gla) residues (Price 
& Williamson, 1985; Loeser & Wallin, 1992) . The gla 
residues are presumably products of a post- translational 
modification by the vitamin K dependent enzyme y- 
carboxylase . MGP strongly binds to hydroxyapatite in a gla 

15 dependent fashion (Dowd et al . , 1995). High levels of MGP 
are found in the extracellular matrix of bone, dentin and 
cartilage (Hale et al . , 1988) . However, MGP is expressed in 
many tissues, with the highest levels of mRNA found in 
lung, heart, kidney and cartilage (Fraser & Price., 1988) . 

2 0 A first indication of the function of MGP in bone came from 
experiments with rats treated with the y-carboxylase 
inhibitor warfarin. These animals showed excessive 
mineralization of the growth plate, indicating that one 
function of MGP in bone and cartilage could be inhibition 

25 of hydroxyapatite formation (Price et al . , 1982). Final 
proof for this function has come from a recent study on MGP 
knockout mice, which die within 2 months after birth as a 
result of arterial calcification leading to blood-vessel 
rupture. MGP deficient mice also display inappropriate 

30 calcification of the growth plate cartilage, where 
calcification has extended into the zone of proliferating 
chondrocytes rather than being restricted to the lower 
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hypertrophic zone as observed in normal animals. The 
abnormal calcification led to a disorganization of 
chondrocyte columns, eventually resulting in short stature, 
osteopenia and fractures (Luo et al . , 1997), These results 
5 strongly suggest that MGP functions to inhibit 
calcification in soft tissues and restrict mineralization 
within the growth plate cartilage to the lower hypertrophic 
zone - the latter possibly by inhibiting calcification in 
the underlying area of proliferating chondrocytes. 

10 Osteopontin is a phosphorylated and glycosylated 

protein of 44 kDa (Prince et al . , 1987). Together with BSP, 
OPN is the most abundant non-collagenous protein in the 
bone matrix (Nagata et al . , 1991), but it is, unlike BSP, 
also expressed in several other tissues (Denhardt & Guo 

15 1993) . Osteopontin and BSP are clearly related: 1) they 
both have an RGD domain that mediates cell attachment via 
ce^&2 integrin class of cell surface receptor (Ross et al . , 
1993; Wong et al . , 1996), and 2) they both have a high 
content of acidic amino acids (OPN has poly-aspart ic acid 

2 0 segments and BSP has poly-glutamic acid segments) and 
sialic acid (Franzen & Heinegard) . 

Phosphorylated osteopontin is a potent inhibitor of 
hydroxyapatite formation while the dephosphorylated form is 
far less potent (Hunter et al . , 1994). In bone osteopontin 

25 is found at high concentrations in the lamina limitans that 
underlies bone lining cells and in reversal (cement) lines 
found at matrix-matrix interfaces where bone deposition has 
been preceded by a resorptive event (McKee & Nanci 1996; 
McKee et al . , 1993) . These findings together with the 

30 hydroxyapatite inhibiting activity of osteopontin suggest 
that osteopontin may act to seal off growing hydroxyapatite 
surfaces once active bone formation has ceased, as 
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speculated by Hunter et al . , 1994. Recently, it has been 
demonstrated that osteopontin knockout mice show normal 
development and bone structure but osteoclast formation is 
enhanced in vitro (Rittling et al . , 1998) and osteoclast 
5 numbers are higher in epiphyseal regions in OPN -/- than in 
wild type mice (Yoshitake et al . , 1998). Interestingly, 
Yoshitake and co-workers also showed that newly formed 
excessive bone following bone marrow ablation in femur from 
wild type animals was resorbed after 2 weeks while no bone 

10 resorbtion was observed in a similar experiment in OPN -/- 
mice (Yoshitake et al . , 1998). Hence, the increased number 
of osteoclasts in OPN -/- mice is likely a compensation for 
the reduced resorbing ability. This substantiates the 
notion put forth by Hunter and co-workers above, that 

15 osteopontin acts to limit growing hydroxyapatite surfaces 
following bone formation. 

Another condition, which has been suggested to be 
promoted by matrix gla protein and osteopontin is 
atherosclerosis {Shanahan et al . , 1994; Sohma et al . , 

20 1994) . Employing differential hybridization techniques to 
screen a cDNA library derived from the aortae of Watanabe 
heritable hyperlipidemic (WHHL) rabbits, Sohma and co- 
workers found one clone encoding the matrix gla protein 
(Sohma et al . , 1994). Northern blot analysis of RNA 

25 prepared from aortae of WHHL and normal rabbits of various 
ages indicated that the expression of matrix gla protein 
mRNA increased in proportion to the progression to 
atherosclerosis in the WHHL rabbits (Sohma et al . , 1994). 
Additionally, mice lacking matrix gla protein die from 

30 aortic calcification (Luo et al . , 1997). No atherosclerotic 
plaques were found in these mice, suggesting that matrix 
gla protein may instead affect calcification once plaques 
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are formed. Nevertheless, high levels of matrix gla protein 
have been located in lipid-rich areas of atherosclerotic 
plaques (Shanahan et al . , 1994). Also, high levels of 
osteopontin ttiRNA and protein have been found in necrotic 
5 lipid cores and areas of calcification in human 
atheromatous plaques (Shanahan et al . , 1994). 

Osteoprotegerin (OPG) /osteolastogenesis inhibitory 
factor (OCIF) was recently identified independently by two 
groups as a 3 80 amino acid residue long glycoprotein with a 

10 molecular weight of approximately 55 kD related to the 
tumor necrosis factor receptor (TNF-R) superfamily (Simonet 
et al., 1997: Yasuda et al . , 1997). Unlike the other TNF- 
R-like molecules this cytokine receptor lacks a trans- 
membrane domain (Simonet et al . , 1997). Accordingly, 

15 OPG/OCIF is a secreted protein appearing as a disulfide 
linked homodimer with a molecular weight of approximately 
110 kD (Simonet et al . 1997). 

Initially, mice transgenic for rat OPG/OCIF were found 
to develop osteopetrosis, which indicated that OPG/OCIF 

2 0 either could function to increase osteoblast -mediated bone 

formation or to decrease osteoclast -mediated bone 
resorption (Simonet et al . , 1997). In an in vitro 

osteoclast -forming assay recombinant OPG/OCIF was found to 
be a potent inhibitor of osteoclastogenesis (Simonet et 
25 al., 1997). In agreement with this finding OPG/OCIF 
knockout mice develop osteoporosis, emphasising that 
OPG/OFIC is an important regulator of postnatal bone mass 
(Bucay et al , , 1998) . 

During mouse embryogenesis OPG/OCIF is highly 

3 0 expressed in cartilage of developing bones, as well as in 

several major arteries, the gastrointestinal tract, and 
skin (Simonet et al . , 1997). In the adult animal OPG/OCIF 
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expression is found in several tissues including heart, 
brain, lung and liver, which is in contrast to human tissue 
where OPG/OCIF expression is absent in brain and liver, 
highly expressed in kidney and detectable in various 
5 hematopoietic and immune organs (Simonet et al . , 1997). No 
explanations have been offered on these observations, 
except that it could be due to true species-specific 
expression differences (Simonet et al . , 1997). 

The ligand for OPG/OCIF has also been identified 

10 independently by the same two groups, who identified 
OPG/OCIF. The ligand, called OPG ligand (OPGL) (Lacey et 
al., 1998) or osteoclast differentiation factor (ODF) 
(Yasuda et al . , 1998), is a TNF-related cytokine, which 
binds to a hematopoietic progenitor cell committed to the 

15 osteoclast lineage and stimulates its differentiation into 
an osteoclast (Lacey et al . , 1998). In vitro OPGL/ODF also 
stimulates mature osteoclasts to resorb bone and 
recombinant OPGL/ODF injected subcutaneously stimulates 
bone resorption in mice (Lacey et al . , 1998). OPGL/ODF is 

20 produced either as a 45 kDa membrane bound protein or as a 
31 kDa soluble, secreted C-terminal fragment (Lacey et al . , 
1998) . 

OPGL/ODF, which is identical to two previously 
identified cytokines, TNF-related activation- induced 

25 cytokine (TRANCE) (Wong et al . , 1997) essential for T-cell 
activation and receptor activator of NF-kB ligand (RANKL) 
essential for dendritic cell activation (Anderson et al . , 
1997) , is highly expressed in lymphoid tissues and 
trabecular bone (Lacey et al . , 1998: Yasuda et al . , 1998). 

30 Since OPG/OCIF binds to and inhibits the action of 

OPGL/ODF, these two proteins are seemingly important extra- 
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cellular regulators of osteoclast development and, thus, 
eventually bone resorption. 

Apart from an osteoporotic phenotype OPG/OCIF knockout 
mice also displayed a marked calcification of the aorta and 
5 renal arteries by 2 months of age (Bucay et al . , 1998) . 
Thus, OPG/OCIF inhibits decalcification of the skeleton and 
at the same time inhibits calcification of certain blood 
vessels. A similar phenomenon has been observed previously 
in matrix gla protein (MGP) knockouts (Luo et al . , 1997). 

10 However, the arterial calcification is more disseminated 
and pronounced in the MGP knockouts, while the bone loss is 
more severe in the OPG/OCIF knockouts. Aberrant expression 
of proteins involved in the prevention of the calcification 
of blood vessels could, theoretically, be associated with 

15 the generation of atherosclerotic plaques. However, no 
atherosclerotic plaques were found in the OPG/OCIF knockout 
mice (Bucay et al . , 1998), excluding a direct role in the 
generation of atherosclerosis. Nevertheless, it would not 
exclude that aberrant expression of OPG/OCIF could be a 

20 secondary factor accelerating the formation of an 
atherosclerotic plaque once the initial • lesion has 
occurred . 

The method of assessing an individual, predisposition 
to osteoporosis or other calcification condition related 

25 diseases described above may be combined with measurements 
of bone mass on a whole body or selected location basis. 
These include X-ray or ultrasound BMD measurements. The 
methods described herein may also or instead be combined 
with measurements of chemical bone resorption markers such 

30 as the Crosslaps™ measurement of C-telopeptide fragments of 
Type 1 collagen or measurements of N-telopeptide fragments 
of Type 1 collagen in body fluids such as serum or urine. 
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Each of these types of measurement may be treated 
as a risk factor to be combined in a weighted manner with 
the one or more of the others (one of course being a 
genetic predisposition measurement according to this 
invention) . 

In order that the nature of the present invention be 
more clearly understood, there now follows an example in 
which reference is made to the Figures shown in the 
accompanying drawings in which: - 

Figure. 1 panel A shows the location of the two 
polymorphisms, called BSP-A1496G and BPS-G1869A, in the 
bone sialoprotein gene promoter. Panel B shows the location 
of the polymorphism, called MGP-C242A, in the matrix gla 
protein gene promoter. Panel C shows the location of the 
5 polymorphisms, called OPN-G520A and OPN-T182 5C, in the 
osteopontin gene promoter. The wild type sequences 
encompassing the four polymorphic sites for all three said 
promoters are shown with the nucleotide at the polymorphic 
position in bold and with the substituting nucleotide - 
0 also in bold - positioned above the polymorphic site. All 
nucleotide numbering is relative to base pair number l, 
which is the most 5' nucleotide of each of the promoter 
sequences as published in the GenBank nucleotide database. 

Figure 2 shows the time dependence of the mean BMCs 
5 measured at the distal arm when grouped into two genotypes 
for the BSP-A1496G polymorphic site. The BMCs were 
determined at four time points - 1977, 1979, 1989, and 1995 
from each individual out of 133 in the 18 years study. Each 
point on the curves is the mean BMC for a given year and a 
0 given genotype with SEM error bars. The table below the 
chart lists the actual BMC values determined in 1977, 1979, 
1989, and 1995 for the two genotype groups as well as 
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difference - in percent - between the two genotype groups 
for each of the four sampling years. 

Figure 3 shows the time dependence of the mean BMCs 
measured at the distal arm when grouped into two genotypes 
5 for the BSP-G1869A polymorphic site. The BMCs were 
determined at four time points - 1977, 1979, 1989, and 1995 
from each individual out of 13 3 in the 18 years study. Each 
point on the curves is the mean BMC for a given year and a 
given genotype with SEM error bars. The table below the 

10 chart lists the actual BMC values determined in 1977, 1979, 
193 9, and 1995 for the two genotype groups as well as 
difference - in percent - between the two genotype groups 
for each of the four sampling years. 

Figure 4 shows time dependence of the mean BMCs 

15 measured at the distal arm when grouped into two genotypes 
for a combination of the BSP-A1496G and BSP-G1869A 
polymorphic sites. The upper curve represents the 
combination: BSP-A14 96G heterozygous /homozygous polymorphic 
and BSP-G1869A wild type while the lower curve represents 

20 the combination: BSP-A1496G wild type and BSP-G1869A 
heterozygous/homozygous polymorphic. The ' BMCs were 
determined at four time points - 1977, 1979, 1989, and 1995 
from each individual out of 133 in the 18 years study. Each 
point on the curves is the mean BMC for a given year and a 

25 given genotype with SEM error bars. The table below the 
chart lists the actual BMC values determined in 1977, 1979, 
1989, and 1995 for the two genotype groups as well as 
difference - in percent - between the two genotype groups 
for each of the four sampling years. 

30 Figure 5 shows the time dependence of the mean BMCs 

measured at the distal arm when grouped into two genotypes 
for the MGP-C242A polymorphic site. The BMCs were 
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determined at four time points - 1977, 1979, 1989, and 1995 
from each individual out of 13 3 in the 18 years study. Each 
point on the curves is the mean BMC for a given year and a 
given genotype with SEM error bars . The table below the 
chart lists the actual BMC values determined in 1977, 1979, 
1989, and 1995 for the two genotype groups as well as 
difference - in percent - between the two genotype groups 
for each of the four sampling years. 

Figure 6 shows the time dependence of the mean BMCs 
measured at. the distal arm when grouped into two genotypes 
for the OPN-G520A polymorphic site. The BMCs were 
determined at four time points - 1977, 1979, 1989, and 1995 
from each individual out of 13 3 in the 18 years study. Each 
point on the curves is the mean BMC for a given year and a 
given genotype with SEM error bars. The table below the 
chart lists the actual BMC values determined in 1977, 1979, 
1989, and 1995 for the two genotype groups as well as 
difference - in percent - between the two genotype groups 
for each of the four sampling years. 

Figure 7 shows the time dependence of the mean BMCs 
measured at the distal arm when grouped into two genotypes 
for a combination of the MGP-C242A and OPN-G520A 
polymorphic sites. The BMCs were determined at four time 
points - 1977, 1979, 1989, and 1995 from each individual 
out of 133 in the 18 years study. Each point on the curves 
is the mean BMC for a given year and a given genotype with 
SEM error bars. The table below the chart lists the actual 
BMC values determined in 1977, 1979, 1989, and 1995 for the 
two genotype groups as well as difference - in percent - 
between the two genotype groups for each of the four 
sampling years. 
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Figure 8 shows the time dependence of the mean BMCs 
measured at the distal arm when grouped into two genotypes 
for the OPN-T1825C polymorphic site. The BMCs were 
determined at four time points - 1977, 1979, 1989, and 1995 
5 from each individual out of 133 in the 18 years study. Each 
point on the curves is the mean BMC for a given year and a 
given genotype with SEM error bars . The table below the 
chart lists the actual BMC values determined in 1977, 1979, 
1989, and 1995 for the two genotype groups as well as 

10 difference - in percent - between the two genotype groups 
for each of the four sampling years. 

Figure 9 shows the time dependence of the mean BMCs 
measured at the distal arm when grouped into two genotypes 
for a combination of the BSP-G1869A and OPN-T1825C 

15 polymorphic sites. The BMCs were determined at four time 
points - 1977, 1979, 1989, and 1995 from each individual 
out of 13 3 in the 18 years study. Each point on the curves 
is the mean BMC for a given year and a given genotype with 
SEM error bars. The table below the chart lists the actual 

20 BMC values determined in 1977, 1979, 1989, and 1995 for the 
two genotype groups as well as difference - ' in percent - 
between the two genotype groups for each of the four 
sampling years. 

Figure 10 shows the location of the OPG-A163G 

25 polymorphism in the OPG/OCIF gene promoter. The wild type 
sequence encompassing the polymorphic site for said 
promoter is shown with the nucleotide at the polymorphic 
position in bold and with the substituting nucleotide - 
also in bold - positioned above the polymoirphic site. All 

30 nucleotide numbering is relative to base pair number 1, 
which is the 5' most nucleotide of the promoter sequence as 
published in the GenBank nucleotide database. 
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Figure 11 shows the time dependence of the mean BMCs 
measured at the distal arm when grouped into two genotypes 
for the OPG-A163G polymorphic site. The BMCs were 

determined at four time points - 1977, 1979, 1989 and 1995 
5 from each out of 133 individuals in the 18 years study. 
Each point on the curves is the mean BMC for a given year 
and a given genotype with SEM error bars . The table below 
the chart lists the actual BMC values determined in 1977, 
1979, 1989 and 1995 for the two genotype groups as well as 

10 difference .- in percent - between the two genotype groups 
for each of the four sampling years. 

Figure 12 shows time dependence of the mean BMCs 
measured at the distal arm when grouped into two genotypes 
for a combination of the BSP-A14 96G and OPG-A163G 

15 polymorphic sites. The upper curve represents the 

combination: BSP-A1496G heterozygous /homozygous polymorphic 
and OPG-A163G wild type while the lower curve represents 
the combination: BSP-A14 96G wild type and OPG-A163G 
heterozygous/homozygous polymorphic . The BMCs were 

20 determined at four time points - 1977, 1979, 1989 and 1995 
from each out of 133 individuals in the IB" year study. 
Each point on the curves is the mean BMC for a given year 
and a given genotype with SEM error bars . The table below 
the chart lists the actual BMC values determined in 1977, 

25 1979, 198 9 and 1995 for the two genotype groups as well as 
difference - in percent - between the two genotype groups 
for each of the four sampling years. 

Figure 13 shows time dependence of the mean BMCs 
measured at the distal arm when grouped into two genotypes 

30 for a combination of the BSP-G1869A and OPG-A163G 
polymorphic sites. The upper curve represents the 

combination: BSP-A1496G homozygous polymorphic and OPG- 
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A163G wild type while the lower curve represents the 
combination: BSP-G1869A wild type/heterozygous and OPG- 
A163G heterozygous/homozygous polymorphic. The BMCs were 
determined at four time points - 1977, 1979, 1989 and 1995 
5 from each out of 133 individuals in the 18 year study. 
Each point on the curves is the mean BMC for a given year 
and a given genotype with SEM error bars. The table below 
the chart lists the actual BMC values deterroined in 1977, 
1979, 1989 and 1995 for the two genotype groups as well as 
10 difference - in percent - between the two genotype groups 
for each of the four sampling years. 

Example 1 (18 year study) 

Methods 

15 Subjects. One hundred thirty three women followed up for 18 
years (1977-1995) with respect to BMD, biochemical markers, 
height, and weight were used in the study. A detailed 
description of the cohort has been previously published 
(Jorgensen et al . , 1996). 

20 

DNA analyses. Screening for the BSP-A1496G 'and the BSP- 
G1869A polymorphisms (basepair numbering according to 
numbering of BSP promoter sequence submitted to GenBank, 
accession #L24756) , the MGP-C242A polymorphism (basepair 

25 numbering according to numbering of MGP promoter sequence 
submitted to GenBank, accession #M55270) , as well as the 
OPN-G520A and OPN-T182 5C polymorphisms (basepair numbering 
according to numbering of osteopontin promoter sequence 
submitted to GenBank, accession #D14813) were performed as 

30 follows: 

The polymerase chain reaction (PGR) was used to 
amplify approximately 250 bp long DNA fragments of the BSP, 
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MGP, and OPN promoters encompassing the BSP-A1496G, BSP- 
G1869A, MGP-C242A, OPN-G520A, and OPN-T1825C polymorphic 
basepairs, PGR techniques are well known in the art and it 
would be within the ambit of a person of ordinary skill in 
5 this art to identify primers for amplifying a suitable 
section of the BSP, MGP and OPN genes including the 
positions 1496bp and 1869bp in the BSP promoter, the 
position 242bp in the MGP promoter, and the positions 520bp 
and 1825bp in the osteopontin promoter. PGR techniques are 

10 described for example in patents US4683202 or EP0200362B1. 
Two hundred ng of genomic DNA was added to 25 jjl reaction 
containing Ix Taq polymerase buffer with 1.5 mM MgCls 
(Perkin Elmer) , 5 nmol of each dNTP, 2 0 pmol of forward and 
reverse primer, and 1.25 units of AmpliTaq Gold (Perkin 

15 Elmer) . The reaction was heated to 95*^C for 9 minutes 
followed by 35 cycles of 95^0 for 3 0 seconds, 460c (BSP- 
A1496G and BSP-G1869A polymorphisms) or 49^C (MGP-C242A 
polymorphism) or 46°C (OPN-G52 0A polymorphism) or 4 8^0 
(OPN-T182 5C) for 3 0 seconds and 12^C for 3 0 seconds - the 

20 latter incubation with a 5 second time extension per cycle. 
The reaction was finally incubated 7 minutes at 72°C for 
completion of the extension reaction. Primer sequences for 
PGR amplification of DNA fragments encompassing the BSP- 
A1496G, BSP-G1869A, MGP-C242A, and OPN-G520A and OPN-T1825C 

25 polymorphic basepairs were: 

BSP-A1495G polymorphism primer set: 
Forward primer : 5 ' - GAA AAG ATA TAT ATA GAA GCC CAA G - 3 ' 
(SEQ ID No. 1) 

Reverse primer : 5 ' - TAA TAT CAT TTG ATG TTT CCT CCT G - 3 ' 
30 (SEQ ID No. 2) 
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BSP-G1869A polymorphism primer set: 
Forward primer: 5'- TTC TTT CGA CAT AGT GAA AAC ACG T - 3' 
(SEQ ID No. 3) 

Reverse primer: 5'- CGT GGA TTC TCA CCA GAA AAC - 3' (SEQ 
5 ID No. 4) 

MGP-C242A polymorphism primer set: 
Forward primer: 5'- CAG TGA GAA AGC TCA TCA CTT GGT C - 3' 
(SEQ ID No. 5) 

Reverse primer: 5'- ATT CTC CCA TCC ATC CAT CCA TGC A - 3' 
10 (SEQ ID No.. 6) 

OPN-G520A polymorphism primer set: 
Forward primer: 5'- CGC TGG AAT TAA GAA AAT TGG TAG A - 3' 
(SEQ ID No. 7) 

Reverse primer: 5'- GTT GTC AAT TTA GTG GAG GGA GAT C - 3' 
15 (SEQ ID No. 8) 

OPN-T1825C polymorphism primer set: 
Forward primer: 5'- GAG TAG TAA AGG ACA GAG GCG AGC T - 3' 
(SEQ ID No. 9) 

Reverse primer: 5' - CTA GCT TTT TCA TTT ACG GGA TGG G - 3' 
2 0 (SEQ ID No. 10) 

To determine the presence or absence of a polymorphic 
genotype of PGR amplified DNA fragments using the above 
mentioned PGR primer sets, restriction enzyme analyses were 

25 performed as follows: 

DNA fragments PGR amplified using the BSP-A1496G 
polymorphism primer set were restricted with Eco T14 I in a 
2 0 }xl reaction containing: Ix buffer H (Amersham 
Pharmacia) , 4 units of Eco T14 I (Amersham Pharmacia) and 5 

30 lil of the cycled PGR reaction. The reaction mixture was 
incubated at 21^C for 1 hour. Four iJ.1 6xgel- loading buffer 
(0.25% bromophenol blue, 0.2 5% xylene cyanol FF, 30% 
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glycerol in water) were added to the 20 ^1 Eco T14 I digest 
and loaded on a 2.5% agarose gel. DNA fragments were then 
resolved by electrophoresis until the bromophenol blue 
marker had run 2/3 through the gel. If the DNA sample 
5 analyzed was homozygous for the wild type BSP sequence one 
band of 270 bp would be observed. If the DNA sample 
analyzed was heterozygous two bands of 270 bp and 245 bp 
would be observed. If the DNA sample was homozygous for the 
polymorphism one band of 245 bp would be observed. 
10 DNA fragments PGR amplified using the BSP-G1869A 

polymorphism primer set were restricted with Eco 72 I in a 
2 0 ^1 reaction containing: Ix Universal buffer 
(Stratagene) , 4 units of Eco 72 I (Stratagene) and 5 /xl of 
the cycled PGR reaction. The reaction mixture was incubated 
15 at 37OC for 1 hour. Four ^1 6xgel - loading buffer (0.25% 
bromophenol blue, 0.25% xylene cyanol FF, 30% glycerol in 
water) were added to the 20 m1 Eco 72 I digest and loaded 
on a 2.5% agarose gel. DNA fragments were then resolved by 
electrophoresis until the bromophenol blue marker had run 
20 2/3 through the gel. If the DNA sample analyzed was 
homozygous for the wild type BSP sequence one band of 2 53 
bp would be observed. If the DNA sample analyzed was 
heterozygous two bands of 253 bp and 23 0 bp would be 
observed. If the DNA sample was homozygous for the 
25 polymorphism one band of 23 0 bp would be observed. 

DNA fragments PGR amplified using the MGP-C242A 
polymorphism primer set were restricted with Eco T22 I in a 
2 0 ^1 reaction containing: Ix buffer H (Amersham 
Pharmacia) , 4 units of Eco T22 I (Amersham Pharmacia) and 5 
30 ^1 of the cycled PGR reaction. The reaction mixture was 
incubated at 37^C for 1 hour. Four ptl 6xgel-loading buffer 
(0.25% bromophenol blue, 0.25% xylene cyanol FF, 30% 
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glycerol in water) were added to the 2 0 /il Eco T22 I digest 
and loaded on a 2.5% agarose gel. DNA fragments were then 
resolved by electrophoresis until the bromophenol blue 
marker had run 2/3 through the gel. If the DNA sample 
analyzed was homozygous for the wild type MGP sequence one 
band of 2 66 bp would be observed. If the DNA sample 
analyzed was heterozygous two bands of 2 66 bp and 241 bp 
would be observed. If the DNA sample was homozygous for the 
polymorphism one band of 241 bp would be observed. 

DNA fragments PGR amplified using the OPN-G520A 
polymorphism primer set were restricted with Bgl II in a 2 0 
jxl reaction containing: Ix buffer H {Amersham Pharmacia) , 4 
units of Bgl II {Amersham Pharmacia) and 5 jil of the cycled 
PGR reaction. The reaction mixture was incubated at 37°C 
for 1 hour- Four ^1 6xgel -loading buffer (0.25% bromophenol 
blue, 0.25% xylene cyanol FF, 3 0% glycerol in water) were 
added to the 2 0 /xl Bgl II digest and loaded on a 2.5% 
agarose gel. DNA fragments were then resolved by 
electrophoresis until the bromophenol blue marker had run 
2/3 through the gel. If the DNA sample analyzed was 
homozygous for the wild type OPN sequence one band of 278 
bp would be observed. If the DNA sample analyzed was 
heterozygous two bands of 2 78 bp and 257 bp would be 
observed. If the DNA sample was homozygous for the 
polymorphism one band of 2 57 bp would be observed. 

DNA fragments PGR amplified using the OPN-T1825C 
polymorphism primer set were restricted with Sac I in a 20 
1x1 reaction containing: Ix buffer L (Amersham Pharmacia), 4 
units of Sac I (Amersham Pharmacia) and 5 /xl of the cycled 
PGR reaction. The reaction mixture was incubated at 31°C 
for 1 hour. Four (xl 6xgel -loading buffer (0.25% bromophenol 
blue, 0.25% xylene cyanol FF, 30% glycerol in water) were 
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added to the 20 iil Sac I digest and loaded on a 2.5% 
agarose gel. DNA fragments were then resolved by 
electrophoresis until the bromophenol blue marker had run 
2/3 through the gel. If the DNA sample analyzed was 
5 homozygous for the wild type OPN sequence one band of 2 56 
bp would be observed. If the DNA sample analyzed was 
heterozygous two bands of 256 bp and 235 bp would be 
observed. If the DNA sample was homozygous for the 
polymorphism one band of 235 bp would be observed. 

10 

Statistical Methods. To test whether the difference between 
mean BHDs of two genotypes was statistically significant 
two-way, unpaired t- tests were applied. 

15 Results 

Five previously unknown polymorphisms were identified 
by sequencing specific promoter regions from the human BSP 
gene promoter, the human MGP gene promoter, and the human 

2 0 OPN gene promoter following a PCR amplification of 40 DNA 
samples from healthy women. The BSP-A1496G, BSP-G18 69A, 
MGP-C242A, OPN-G520A, and OPN-T1825C polymorphisms were 
coded as Xx, Yy, Zz, Bb, and Ss, respectively, where the 
uppercase letter signifies presence of the wild type base 

25 pair at the given polymorphic position and the lowercase 
letter signifies presence of the base pair different from 
the wild type base pair at the given polymorphic position. 

The 133 DNA samples from the 18 year study were 
screened for the presence of either of the 5 polymorphisms. 

30 Table 1 shows the genotype distribution of DNA samples 

from the 18 years study for all identified polymorphic 
sites: BSP-A1496G, BSP-G1869A, MGP-C242A, OPN-G520A, and 
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OPN-T1825C. The left panel shows the actual number of 
samples categorized into three genotypes for the 3 
identified polymorphic sites. The right panel displays the 
same analysis as the left except that the numbers represent 
the percent of total samples analyzed for each polymorphic 
site . 

wt = XX, YY, ZZ, BB or SS 
hz = Xx, Yy, Zz, Bb or Ss 
pm = XX, yy, zz, bb or ss 
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18 Year Study 
Genotype distribution 



Number of 
samples 



% of total 



BSP-A1496G 



BSP-G1869A 



MGP-C242A 



OPN-G520A 



0PN-T1825C 



Wt hz pm total wt hz pm Total 

127jBSP-A1496G 



3 BSP-G1869A 



1 MGP-C242A 



3 OPN-G520A 



119 OPN-T1825C 



20.2 100.0 



The genotype distributions for the 5 polymorphisms are 
shown in table 1. For the BSP-A1496G polymorphism the 
homozygous wild type genotype was the most abundant, 
followed by the heterozygous and homozygous polymorphic 
genotypes, with the homozygous polymorphic groups being 
quite small. In the case of the BSP-G1869A polymorphism, 
the wild type genotype, as defined by the BSP gene promoter 
sequence from GenBank, was rare, the heterozygous genotype 
was 10 times more frequent and the homozygous polymorphic 
genotype was twice as frequent as the heterozygous. For the 
MGP-C242A polymorphism the heterozygous genotype was the 
most abundant followed by the homozygous wild type and 
homozygous polymorphic genotypes . In the case of the OPN- 
G52 0A polymorphism, the homozygous polymorphic genotype was 
the most abundant followed by the heterozygous and the wild 
type homozygous genotypes. The genotype distribution of the 
OPN-T182 5C polymorphism was, generally, the same as for the 
MGP-C242A polymorphism. 
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The impact of the 5 identified polymorphic sites on 
bone mass as represented by BMC and BMD measurements at the 
distal arm, as well as percent change of BMC/BMD over time, 
was analyzed. A compilation of these analyses is shown in 
5 table 2 . 

Table 2 shows a compilation of a statistical analysis 
of the results obtained from screening of the DNA samples 
from the 18 years study for the presence of either of the 5 
identified polymorphisms. Panel A: The numbers represent 

10 the likelihood that the difference in mean BMD between 
groups of different genotype {homozygous wild type or 
heterozygous/homozygous polymorphic) are identical. The 
test includes BMC {bone mineral content) and BMD measured 
at the distal arm. Numbers in parentheses represent the 

15 year of BMC/BMD measurement- Panel B: The numbers represent 
the difference between the mean HMDs of two genotype groups 
{heterozygous/homozygous polymorphic group subtracted from 
the homozygous wild type group) in percent of the highest 
BMD for a given polymorphic site. 
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From this table it is clear that the BSP-A1496G and 
BSP-G1869A polymorphic sites, especially when combined, are 
good sites for predicting whether an individual is 
genetically predisposed for high or low BMC/BMD. The OPN- 
T1825C polymorphism only has a marginal influence on 
BMC/BMD on its own. However, when the OPN-T182 5C 
polymorphism is combined with the BSP-G1869A polymorphism 
the percent separation of genotypes is better than either 
polymorphism alone (Table 2) . On the other hand, the MGP- 
C242A and OPN-G520A polymorphisms are, at first glance, not 
suitable sites for such a prediction. None of the 
identified polymorphisms appeared to have a statistically 
significant impact on the change in bone mass over time 
{data not shown) . Age, height and weight of the individuals 
involved in the 18 years study did not differ significantly 
between any of the genotype groups (data not shown) . 

These observations strongly indicated that the ESP 
polymorphisms influence peak bone mass rather than the rate 
of bone loss. To substantiate this an analysis of the 
variation of BMD as measured 4 times on the same individual 
from 1977 to 1995 for the different genotypes was 
performed. In 1977 the average age of the individuals 
included in the 18 years study was 51.1 years, thus ending 
at 69.1 years in 1995. The expected outcome of a plot of 
the means of BMD for one of the two BSP polymorphisms as a 
function of time would be two parallel curves, each 
representing HMDs measured in individuals with the wild 
type genotype and BMDs measured in individuals with the 
polymorphic phenotype. Figures 2 and 3 show that this is, 
indeed, the case for the BSP-A1496G and BSP-G1869A 
polymorphic sites. Moreover, the two BSP promoter 
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polymorphisms act in concert on peak bone mass to augment 
the mean BMD difference between genotypes even more than 
the isolated contribution of each polymorphism (Figure 4) . 

The role - if any - of the MGP-C242A and OPN-G520A 
5 polymorphisms in the MGP and OPN promoters on bone turnover 
was less clear from the first analyses compiled in table 2 . 
However, when BMC values grouped according to genotype were 
plotted as a function of time a set of curves appeared 
suggesting that both the MGP-C242A (Figure 5) and OPN-G52 0A 

10 (Figure 6) polymorphic sites are determinants of rate of 
bone loss. It is especially noteworthy that the ZZ and 
Zz+22 curves as well as the BB+Bb and bb curves separate 
between 1979 and 1989, corresponding to an average age of 
53.1 years and 63,1 years, indicative of a genetic 

15 phenomenon associated with the menopause. Like the BSP 
polymorphisms, the combined action of the MGP-C242A and 
OPN-G52 0A polymorphisms also leads to a bigger difference 
between genotypes than either would create alone (Figure 
7) . 

20 According to the results compiled in table 2, the 

impact of the OPN-T1825C polymorphism on BMC/BMD was only 
visible after it was combined with the BSP-G1869A 
polymorphism. From a graph of BMC values grouped according 
to genotype and plotted as a function of time it is 

25 difficult to tell whether this polymorphism has an impact 
on rate of bone loss or peak bone mass, due to the 
proximity of the curves (Figure 8} . However, combining this 
polymorphism with the BSP-G18 69A polymorphism gave rise to 
a set of time course curves clearly showing that these two 

3 0 polymorphisms cooperate in an additive fashion, and, hence, 
that the OPN-T182 5C polymorphism may influence peak bone 
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mass, as the BSP-G1869A polymorphism, rather than the rate 
of bone loss (Figure 9) . 

Finally, the association between genotype and urinary 
osteocalcin (N-MID®, Osteometer Biotech A/S) as well as 
5 urinary collagen type 1 C-terminal crosslinks (CrossLaps®, 
Osteometer Biotech A/S) were examined. As expected there 
was no significant difference between the mean value of 
either of the biochemical bone turnover markers for the 
BSP-A1496G and BSP-G1869A polymorphic sites (data not 
10 shown) . Also, no significant difference was observed 
between the mean value of the biochemical bone turnover 
markers for the MGP-C242A and OPN-G520A polymorphic sites. 
This is likely due to the equally paced bone loss of the ZZ 
^ and Zz+zz genotype groups as well as the BB+Bb and bb 

^ 15 genotypes groups at the time of N-MID® and CrossLaps® 

measurement (1995) according to Figures 5 and 6. 

Example 2 (18 year study) 

Methods 

20 Subjects. One hundred thirty three women followed up for 
18 years (1977-1995) with respect to bone mineral content 
(BCM) or bone mineral density (BMD) , biochemical markers, 
height, and weight were used in the study. A detailed 
description of the cohort has been previously published 

25 (J0rgensen et al , , 1996). 

DNA analyses. Screening for the OPG-A163G polymorphism 
{basepair numbering according to numbering of OPG/OCIF 
promoter sequence submitted to GenBank, accession 
30 #AB008821) , was performed as follows: 
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The polymerase chain reaction (PCR) was used to 
amplify a 253 bp long DNA fragment of the OPG/OCIF promoter 
encompassing the OPG-A163G polymorphic basepair. Two 
5 hundred ng of genomic DNA was added to 25 jul reaction 
containing IxPCR Gold buffer {Perkin Elmer), 1.5 mM MgCl2,5 
nmol of each dNTP, 2 0 pmol of forward and reverse primer, 
and 1.25 units of AmpliTaq Gold (Perkin Elmer) . The 
reaction was heated to 95°C for 9 minutes followed by 35 
10 cycles of 95°C for 30 seconds, 46°C for 30 seconds and 72°C 
for 3 0 seconds - the latter incubation with a 5 second time 
extension per cycle. The reaction was finally incubated 7 
minutes at 12°C for completion of the extension reaction. 
The primer sequences for this PCR amplification were: 
15 OPG-A163G polymorphism primer set: 

Forward primer: 5 ' -AGT CTA ACT TCT AGA CCA GGC AAT T-3' 
(SEQ ID No. 11) 

Reverse primer: 3 ' -AGT TAG AGC CAG AGA GAA TCT G-3' (SEQ 
ID No. 12) 

20 

To determine the presence or absence of 'a polymorphic 
genotype of PCR amplified DNA fragments using the PCR 
primer sets above, restriction enzyme analyses were 
performed as follows: 

25 DNA fragments PCR amplified using the OPG-C163A 

polymorphism primer set were restricted with Mfe I in a 20 
fil reaction containing: 1 x NEBuffer 4 (New England 
Biolabs) , 4 units of Mfel (New England Biolabs) and 5 fil of 
the cycled PCR reaction. The reaction mixture was 

3 0 incubated at 3 7°C for 1 hour. Four /il 6xgel- loading buffer 
(0.25% bromophenol blue, 0.25% xylene cyanol FF, 3 0% 
glycerol in water) were added to the 20 ^1 Mfe I digest and 
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loaded on a 2.5% agarose gel. DNA fragments were then 
resolved by electrophoresis until the bromophenol blue 
marker had run 2/3 through the gel. If the DNA sample 
analysed was homozygous for the wild type OPG/OCIF sequence 
5 one band of 253 bp would be observed. If the DNA sample 
analysed was heterozygous two bands of 2 53 bp and 23 2 bp 
would be observed. If the DNA sample was homozygous for 
the polymorphism one band of 232 bp would be observed. 

10 Statistical Methods. To test whether the difference 
between mean BMDs of two genotypes was statistically 
significant two-way, unpaired t-tests were applied. 

Results 

15 

A previously unknown polymorphism, OPG-A163G, was 
identified by sequencing the promoter region from the human 
OPG/OCIF gene promoter, following a PGR amplification of 40 
DNA samples from healthy women. The OPG-A163G polymorphism 

2 0 was coded as Mm, where the uppercase letter signifies 
presence of the wild type base pair at the given 
polymorphic position and the lowercase letter signifies 
presence of the base pair different from the wild type base 
pair at the given polymorphic position. The position of 

25 the polymorphism is depicted in Figure 10. 

The 133 DNA samples from the 18 year study were 
screened for the presence of the OPG-A163G polymorphism. 
The genotype distribution of DNA samples from the 18 year 
study was as follows: Homozygous wild type (MM) = 71.3% 

30 (n=92) , heterozygous (Mm) =25.6% (n = 33), and homozygous 
polymorphic (mm) =3.1% (n = 4) . The genotype could not be 
determined in 4 out of the 133 DNA samples. 
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The impact of the identified polymorphic site on bone 
mass as represented by bone mineral content (BMC) and bone 
mineral density (BMD) measurements at the distal arm in 
1977 and 1995, respectively, was analysed (Table 3) . The 
5 percent difference between the genotype groups did not 
change significantly from 1977 to 1995, which implied that 
the OPG-A163G polymorphism exerted an influence on peak 
bone mass. The two polymorphisms, called BSP-A14 96G and 
BSP-G18 69A, described above also have an impact on peak 

10 bone mass." Hence, it was of interest to examine whether 
any co-operation between the OPG polymorphism and either of 
the ESP polymorphisms existed. The combinations OPG- 
A163G/BSP-A1496G and OPG-A163G/BSP-G1869A showed that those 
polymor-phisms certainly act in a co-operative fashion, in 

15 that the t-test p-values for the null -hypothesis (i.e. no 
difference between the genotype groups) dropped to 
statistically significant values and the percent difference 
in mean BMC/BMD values for two genotype groups increased 
(Table 3) . 
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A) 

5 t-test, p- values for the "wild type genotype group = 
heterozygous + homozygous genotype polymorphic groups" 
hypothesis . 





OPG-A163G 


OPG-A163G/ 
BSP-A1496G 


OPG-A163G/ 
BSP-G1869A 


BMC, arm (77) 


0.022 


0.003 


0.056 


BMD, arm (95) 


0.281 


0 . 018 


0 .030 


B) 

BMC and BMD percent differences between genotype groups . 




OPG-A163G 


OPG-A163G/ 
BSP-A1496G 


OPG-A163G/ 
BSP-G1869A 


BMC, arm (77) 


6.7 


12 .3 


8.7 


BMD, arm (95) 


4 . 2 


13 .3 


12 .3 



Thus, it is clear that the OPG-A163G -polymorphism, 
especially in combination with the BSP-A1496G and BSP- 
G1B69A polymorphisms, is a good site for predicting whether 

15 an individual is genetically predisposed for high or low 
BMC/BMD. Age, height and weight of the individuals 
involved in the 18 year study did not differ significantly 
between any of the genotype groups (data not shown) . 

To substantiate the initial indication, that the 

2 0 OPG/OCIF polymorphism influences peak bone mass, an 
analysis of the variation of BMC as measured 4 times on the 
same individual from 1977 to 1995 for the different 
genotypes was performed. In 1977 the average age of the 
individuals included in the 18 year study was 51.1 years, 
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10 



15 



thus ending at 69.1 years in 1995. The expected outcome of 
a plot of the BMC means for the OPG-A163G polymorphism as a 
function of time would be two parallel curves, each 
representing BMCs measured in individuals with the wild 
type genotype and BMCs measured in individuals with the 
heterozygous or polymorphic homozygous phenotypes . Figure 
11 shows that this is certainly the case. Moreover, the 
combination of the OPG-A163G and the BSP-A1496G 
polymorphisms show that they act in concert on peak bone 
mass to augment the mean BMD difference between genotypes 
even more than the isolated contribution of each 
polymorphism (Figure 12). In fact, this co-operation is 
completely additive (Table 4) , indicating that the two 
polymorphisms act on bone mass independently of one 
another. The numbers for the BSP-A1496G and BSP-G1869A 
polymorphisms in Table 4 are from the results presented in 
Figures 2 and 3. Also, the combination of the OPG-A163G 
and BSP-G1869A polymorphisms indicates a positive co- 
operation (Figure 13) , which is almost additive (Table 4) . 



20 



Table 4 



% difference between the "high BMD" and "low BMD" genotype 



groups . 





1977 


1979 


1989 


1995 


BSP-A1496G 


5.1 


8.3 


7.0 


8.5 


BSP-G1869A 


5.0 


5.2 


8 . 8 


8.3 


OPG-A163G 


6.7 


9.2 


4.6 


4.9 


OPG-A163G/BSP-A1496G 


12 .3 


17.8 


12 .7 


13 .5 


OPG-A163G/BSP-G1869G 


8.7 


12 .8 


13 .0 


12.7 
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Finally, the association between genotype and urinary 
osteocalcin (N-MID®, Ostometer Biotech A/S} as well as 
urinary collagen type 1 C-terminal crosslinks (CrossLaps®, 
Osteometer Biotech A/S) were examined. The urine samples 
5 were all collected post the age of attainment of peak bone 
mass. No significant difference between the mean value of 
either of the biochemical bone turnover markers for the 
OPG-A163G polymorphism was found (data not shown) . This 
result is not surprising for a polymorphism with an effect 
10 on peak bone mass, due to the equally paced bone loss of 
the MM and Mm+mm genotype groups at the time of N-MID® and 
CrossLaps® measurement (1995) . 
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A method of assessing an individual's predisposition 
to a selected calcification condition status, which 
method comprises determining the genotype of the 
promoter of the bone sialoprotein gene, the promoter 
of the matrix gla protein gene, the promoter of the 
osteopontin gene, or the osteoprotegerin gene, also 
known as the osteoclastogenesis inhibitory factor gene 
{OPG/OCIF) , or all four or any combination of two or 
more out of the four promoters. 

A method as claimed in Claim 1, wherein said 
calcification condition status is having a high or low 
peak bone mass or having a high or low rate of bone 
loss . 

A method as claimed in Claim 1 or Claim 2, wherein it 
is determined whether the individual is homozygous or 
heterozygous for an allelic variation of the promoter 
of the bone sialoprotein gene, the promoter of the 
matrix gla protein gene, the promoter of the 
osteopontin gene or the promoter of the OPG/OCIF gene, 
or all four or a combination of two or more out of the 
four promoters. 

A method as claimed in Claim 3, wherein said allelic 
variation of the bone sialoprotein gene promoter is 
BSP-A1496G or BSP-G1869A. 

A method as claimed in Claim 3, wherein said allelic 
variation of the matrix gla protein gene promoter is 
MGP-C242A. 

A method as claimed in Claim 3, wherein said allelic 
variation of the osteopontin gene promoter is OPN- 
G52 0A or OPN-T1825C. 

A method as claimed in Claim 3, wherein said allelic 
variation of the osteoprotegerin gene, also known as 
the osteoclastogenesis inhibitory factor gene is OPG- 
A163G. 

A method as claimed in Claim 7, further comprising 
determining whether the individual is homozygous or 
heterozygous for an allelic variation of the promoter 
of the bone sialoprotein gene. 
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9. A method as claimed in Claim 8, wherein said allelic 
variation of the bone sialoprotein gene promoter is 
BSP-A1496G or BSP-G1869A. 

10. A method as claimed in Claim 3, comprising determining 
whether at least one copy of the bone sialoprotein 
gene promoter of the individual has adenine or guanine 
at position 1496 bp or at base 1869 bp, wherein 
adenine at position 1496 bp and guanine at position 
1869 bp are associated with a lower peak bone mass. 

11. A method as claimed in Claim 3, comprising determining 
whether at least one copy of the matrix gla protein 
gene promoter of the individual has cytosine or 
adenine at position 242 bp, wherein adenine is 
associated with a higher rate of loss of bone mass. 

12. A method as claimed in Claim 3, comprising determining 
whether at least one copy of the osteopontin gene 
promoter of the individual has guanine or adenine at 
position 520 bp or thymine or cytosine at position 
1825 bp, wherein adenine at position 520 bp is 
associated with a higher rate of loss of bone mass and 
thymine at position 1825 bp is associated with a lower 
bone mass . 

13. A method as claimed in Claim 3, comprising determining 
whether at least one copy of the osteoprotegerin/ 
osteoclastogenesis inhibitory factor gene promoter of 
the individual has adenine or guanine at position 163 
bp, wherein guanine at position 163 bp is associated 
with a lower peak bone mass . 

14. A method as claimed in any preceding claim, comprising 
amplifying a relevant portion of the DNA of a said 
gene promoter of said individual . 

15. A method as claimed in Claim 14, wherein the sequence 
of said amplified portion is determined by 
hybridisation assay or by restriction fragment length 
analysis . 

16. An oligonucleotide primer for use in amplification of 
a relevant portion of a said gene promoter. 
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DNA comprising a bone sialoprotein gene, or fragment 
thereof at least 15 nucleotides in length, in which 
adenine at position 14 96 bp is substituted by guanine, 
or DNA complementary thereto. 

DNA comprising a bone sialoprotein gene, or fragment 
thereof at least 15 nucleotides in length, in which 
guanine at position 1869 bp is substituted by adenine, 
or DNA complementary thereto. 

DNA comprising a matrix gla protein gene, or fragment 
thereof at least 15 nucleotides in length, in which 
cytosine at position 242 bp is substituted by adenine, 
or DNA complementary thereto. 

DNA comprising an osteopontin gene, or fragment 
thereof at least 15 nucleotides in length, in which 
guanine at position 52 0 bp is substituted by adenine, 
or DNA complementary thereto. 

DNA comprising an osteopontin gene, or fragment 
thereof at least 15 nucleotides in length, in which 
thymine at position 1825 bp is substituted by 
cytosine, or DNA complementary thereto. 

DNA comprising an osteoprotegerin/osteoclastogenesis 
inhibitory factor gene, or fragment thereof at least 
15 nucleotides in length, in which adenine at position 
163 bp is substituted by guanine, or DNA complementary 
thereto . 

A method of osteoporosis therapy comprising 
determining a predisposition as claimed in any one of 
Claims 1 to 22, and administering a medicament to the 
individual to prevent or treat osteoporosis or to 
delay the onset of osteoporosis if the individual is 
predisposed to low peak bone mass or to a high rate of 
loss of bone mass. 

A method of atherosclerosis therapy comprising 
determining a predisposition as claimed in any one of 
Claims 1 to 22, and administering a medicament to the 
individual to prevent or treat atherosclerosis or to 
delay the onset of atherosclerosis if the individual 
is predisposed to pathological arterial calcification. 
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FIGURE 1 
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BSP-A1496G, 18 Years Study 
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FIGURE 3 
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FIGURE 4 





BSP-A1496G/BSP-G1869A, 18 Years Study 
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FIGURE 5 
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